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Retrofit  enclosure  of  gas  furnaces  to  supply  outdoor 
combustion  and  dilution  air  may  save  up  to  6  percent  of 
natural  gas  usage  in  cold  climates-  Retrofit  must  be  in 
accordance  with  local  building  codes  and  the  National  Fuel 
Gas  Code.  NFPA  54.  Cost  effectiveness  of  this  type  of 
energy  conserving  retrofit  will  be  quite  variable  and 
should  be  carefully  calculated  for  each  candidate  building 
before  a  decision  is  reached.  Comparisons  with  other 
retrofit  alternatives  such  as  vent  dampers  and  derating 
should  be  made  before  a  final  retrofit  decision  is  made. 

Gas  water  heaters  in  cold  climates  may  preclude  furnace 
combustion  air  retrofits  because  the  furnace  enclosure 
which  should  include  the  gas  water  heater  will  nearly  reach 
outside  air  temperature  which  could  freeze  pipes  and  cause 
expensive  damage. 

A  natural  gas  furnace  test  facility  at  Colorado  State 
University  was  modified  to  test  the  effect  the  delivery  of 
unheated  outdoor  combustion  and  dilution  air  has  on  thermal 
efficiency.  Tests  were  scheduled  over  the  winter  months  so 
that  the  resulting  efficiency  with  cold  temperatures  in  the 
furnace  room  could  be  determined-  Resulting  efficiencies 
were  compared  to  a  "baseline"  case  where  combustion  and 
dilution  air  were  tempered  indoor  air- 

A  furnace  "model"  produced  from  the  actual  data  was 
interfaced  with  a  computer  simulation  program  called 
"TRNSYS"  to  test  the  effect  of  outdoor  combustion  air  on 
annual  fuel  consumption-  This  "TRNSYS"  simulation  program 
was  developed  by  the  University  of  Wisconsin  at  Madison  and 
uses  actual  weather  data  tapes,  called  "Typical 
Meteorlogical  Year"  data  to  drive  the  models. 


The  data  obtained  from  the  furnace  testing  showed  an 
overall  reduction  in  furnace  cyclic,  thermal  efficiency  as 
outdoor  air  temperature  decreased.  This  reduction  in 
efficiency  was  most  extreme  for  the  short  10  percent  duty 
cycle  and  least  extreme  for  the  100  percent  (steady  state) 
duty  cycle.  Wind  induced  draft  was  found  to  have  no 
appreciable  effect  on  furnace  efficiencies-  However,  the 
computer  simulation  results  indicated  that  annual  fuel 
savings  of  around  5  percent  can  be  achieved  with  outdoor 
combustion  air  despite  the  reduction  in  furnace  thermal  / 
efficiency  with  cold  outdoor  air  temperatures-  These 
savings  were  accomplished  because  of  the  reduced  air 
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infiltration  into  the  house  normally  required  to  replenish 
combustion  and  dilution  air  drawn  up  the  stack  from  indoor 
air.  Also  wind  induced  draft,  even  when  the  furnace  is 
off.  provides  a  further  penalty  since  it  causes  additional 
air  intrusion  into  the  heated  space. 

It  is  also  reasonable  to  assume  that  during  severe  cold 
weather  the  furnace  is  not  operating  at  the  short  duty 
cycles;  and  therefore,  the  reduction  in  efficiency  is  at  a 
minimum.  Conversely,  the  furnace  is  likely  to  operate  for 
short  duty  cycles  when  the  outside  air  temperature  is 
higher  and  the  reduction  in  efficiency  will  be  less  severe 
even  at  the  10  percent  duty  cycle. 

Temperatures  in  the  furnace  room  were  quite  cold,  below 
freezing  at  times,  and  stayed  cold  even  while  the  furnace 
was  firing  at  100  percent  duty  cycles.  Since  the  natural 
inclination  is  to  place  a  gas  fired  hot  water  heater  in  the 
same  enclosed  space  as  the  furnace,  some  care  and  thought 
must  be  exercised.  Any  additional  plumbing  which 
coincidentally  passes  through  the  furnace  space  must  also 
be  either  protected  or  moved  since  freezing  is  a  definite 
possibility. 


PREFACE 

The  laboratory  study  of  gas  furnace  efficiency  with  outside 
combustion  air  was  conducted  at  Colorado  State  University 
by  Nr.  Thomas  E.  Brisbane  assisted  by 

Ms.  Kathleen  T..  Hancock,  both  of  the  Research  Institute  of 
Colorado  staff. 

Computer  analysis  was  accomplished  on  the  Control  Data 
Cyber  172  computer  at  Colorado  State  University  for  both 
data  reduction  and  simulation. 

Appreciation  is  expressed  to  Ms.  Diane  Riggi.  a  student  at 
CSU,  who  gave  willingly  and  generously  of  her  time  and 
energy,  particularly  during  the  most  adverse  of  weather 
conditions. 

Thanks  and  admiration  to  Stuart  Waterbury.  a  graduate 
student  at  the  Solar  Energy  Application  Lab.  CSU  for  the 
ease  with  which  he  was  able  to  model  and  fit  the  furnace 
date  to  an  existing  and  excellent  simulation  program.  His 
intimate  familiarity  with  the  simulation  yielded  the 
essential  results  of  the  study. 
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INTRODUCTION 


Report  FESA-TS-2072 ,  "Energy  Saving  Devices  for  Gas 
Furnaces",  described  the  performances  of  three  types  of 
devices  designed  to  improve  furnace  seasonal  efficiency. 
Report  FESA-TS-2070 .  "Techniques  for  Control  of  Air 
Infiltration  in  Buildings",  suggested  that  provision  of 
outside  combustion  air  could  provide  insurance  against  the 
risks  of  retrofitting  a  house  "too  tightly".  Outside 
combustion  air  was  claimed  by  some  authorities  to  provide 
another  means  for  improving  furnace  seasonal  efficiency. 
Accordingly,  the  objective  of  this  Task  Order  No.  15  was  to 
determine  the  effect  on  furnace  efficiency  when  outside  air 
is  used  for  combustion  and  dilution  air.  The  impact  of 
outside  furnace  air  on  home  energy  use  and  heating  cost  was 
to  be  determined. 

It  has  also  been  proposed  by  some  that  providing  outside 
combustion  air  to  an  enclosed  furnace  room  can 
substantially  reduce  air  infiltration  into  the  heated 
envelope  whether  that  infiltration  is  intentional  to  allow 
for  adequate  combustion  or  unintentional  through  all  the 
typical  cracks  and  leaks  in  a  house.  The  effect  on  furnace 
efficiencies  at  greatly  reduced  outside  temperatures, 
zero  °F  and  below  is  of  particular  interest  as  this  is 
when  the  benefit  of  reducing  air  infiltration  may  be 
greatest.  The  benefit  of  reduced  air  infiltration  into  the 
heated  envelope  may  not  pay  off  if  furnace  efficiencies  are 
drastically  reduced  by  this  low  temperature  combustion  and 
dilution  air.  Therefore,  just  finding  what  effect  the 
outside  combustion  air  has  on  the  furnace  does  not  answer 
the  entire  question.  The  complete  problem  involves  the 
overall  interrelationship  between  efficiencies  and  air 
infiltration. 

In  order  to  obtain  this  information  computer  simulation  of 
a  furnace/house  model  would  be  required  for  both  the  inside 
and  the  outside  combustion  air.  Further,  this 
furnace/house  should  be  operated  in  several  parts  of  the 
country  to  determine  if  an  indicated  solution  only  works  in 
certain  climates  and  not  in  others.  For  example,  the 
outside  combustion  air  may  pay  off  in  milder  climates  but 
might  not  prove  feasible  in  the  cold  northeast. 

Tests  of  furnace  thermal  efficiency  with  outside  air  were 
conducted  over  a  range  of  outside  air  temperatures  and 
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simulated  wind  conditions-  Tests  were  also  made  for 
thermal  efficiency  with  inside  air  provided  to  the  furnace 
room  for  comparison.  Temperatures,  heat  gains,  and  air 
flow  rates  were  among  the  parameters  measured  to  describe 
furnace  performance.  The  wind  which  creates  a  larger  draft 
effect  up  the  stack  was  simulated  for  two  wind  speeds  by 
using  a  two-speed  induced  draft  fan  on  the  stack.  The 
furnace  was  operated  through  a  full  range  of  duty  cycles 
for  each  selected  outdoor  temperature  and  wind  condition. 

Once  the  effect  of  outside  combustion  air  was  known,  the 
furnace  performance  was  simulated  in  a  computer  model  which 
included  a  typical  house  and  appropriate  losses.  The 
computer  models  compared  the  annual  fuel  consumption  for 
two  similar  houses,  one  with  outside  combustion  air  and  one 
with  inside  combustion  air  for  three  locations  around  the 
country.  Actual  meteorological  data  tapes  are  played 
against  the  model  to  determine  this  annual  fuel  utilization 
figure  for  Boston.  MA;  Great  Palls,  MT;  and  Albuquerque. 

NM. 

This  report  describes  the  details  of  the  testing  for 
furnace  efficiency  and  the  results  of  those  tests.  The 
effect  of  outside  combustion  air  on  overall  heating 
performance  is  then  judged  by  comparing  the  annual  fuel 
utilization  using  outside  and  inside  combustion  air  in  a 
typical  house  for  five  different  locations  around  the  U.S. 
This  is  the  result  of  the  output  of  the  computer  simulation 
program  which  models  the  furnace  according  to  the  results 
of  the  earlier  tests.  The  test  facility,  test  procedures 
and  instrumentation  are  described  in  Appendix  A.  Key  data 
are  displayed  in  Appendix  B.  Equations  are  in  Appendix  C 
and  the  computer  program  is  described  in  Appendix  D. 


DISCUSSION 


Testing  Method 

The  test  facility  consisted  of  a  "furnace  room"  6  feet  wide 
and  12  feet  long  and  an  instrumentation  room  the  same  size 
located  adjacent  to  the  furnace  room.  The  "room  air"  for 
the  test  furnace  was  supplied  from  laboratory  space  and 
heated  air  was  delivered  to  the  lab  space  at  8  feet  above 
the  floor  and  directed  in  an  opposite  direction  from  the 
cold  air  intake.  The  air  temperature  near  the  floor  of  the 
laboratory  varied  between  65°F  to  70°F  which  is 
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representative  of  return  air  temperature  in  a  heated 
building. 

A  requlated  natural  gas  supply  line  already  existed  within 
the  building  and  was  easily  extended  to  the  test  rooms  for 
operation  of  the  furnace.  The  combustion  vent  was  run 
through  an  office  balcony  located  immediately  above  the 
test  rooms  and  vented  into  a  space  near  the  "roof"  inside 
the  building.  A  vent  fan.  operating  whenever  the  furnace 
was  operating,  then  pulled  the  combustion  gases  outside  the 
lab  building.  To  simulate  a  wind  induced  draft  through  the 
vent,  a  two-speed  vent  fan  was  used.  Combustion  air  for 
the  furnace  room  was  drawn  through  one  of  two  alternate 
sources:  (1)  for  outside  combustion  air  tests,  the  air  was 

drawn  in  through  holes  drilled  in  the  common  concrete  lab 
wall  directly  to  outside  air-  and  (2)  for  inside  combustion 
air  tests,  tempered  indoor  air  was  drawn  in  from  the  main 
lab  floor. 

The  furnace  was  a  standard  model  which  included  a  four 
speed  blower  and  a  standing  gas  pilot  light  and  all  burner 
controls,  fan  switches  and  24  volt  control  wiring. 

Each  air  duct  monitor  consisted  of  total  and  static 
pressure  probe  arrays  which  were  separately  manifolded  to 
an  outlet  port. 

The  following  data  were  taken  for  the  analysis  of  the 
effect  of  outside  combustion  and  dilution  air  to  the 
furnace  room: 

1-  Temperatures  throughout  the  system- 

2.  Temperature  gains  across  key  elements  of  the  system. 

3.  Air  flow  rates  in  the  room  air  circuit. 

4.  Vent  stack  flow  rates  during  and  before  burner 
operation. 

5.  Wet  and  dry  bulb  temperatures  in  the  laboratory, 
outside  the  furnace  test  room. 

6.  Gas  consumption  during  each  furnace  test. 

7.  Combustion  gas  samples  for  CO2 .  O2  and  CO. 

8.  Natural  gas  samples  for  heating  value  analysis. 
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A  detailed  description  of  testing  apparatus,  the  test 
set-up,  and  data  acquisition  appears  as  Appendix  A. 

The  testing  matrix  for  determining  the  effect  of  outside 
combustion  air  on  furnace  efficiency  was  as  follows: 


Outside  Air 


Temperatures  £F 

Vent  Fan 

Duty  Cycle  Percent 

-10,0,10,25,40,55 

OFF 

10,20,30,50,75,100 

LOW 

10,20,30,50,75,100 

HIGH 

10,20,30,50,75,100 

Indoor  Tempered 

Air 

OFF 

10,20,30,50,75,100 

65  to  70 

LOW 

10,20,30,50,75,100 

HIGH 

10,20,30,50,75,100 

The  low  and  high  stack  fan  settings  were  selected  to 
simulate  the  effect  of  two  different  winds  blowing  on  the 

stack  outlet  on  the  roof  and  inducing  a  draft  in  the 

stack.  The  low  fan  speed  simulated  a  wind  speed  of  about 

25  MPH.  A  high  fan  speed  simulated  a  wind  speed  ranging 

from  40  to  50  MPH.  With  the  stack  fan  off,  the  draft 
effect  was  that  of  the  combustion  process  only  as  if  in 
still  air. 

Data  were  recorded  in  digital  form  on  magnetic  tape.  The 
magnetic  digital  tapes  created  during  each  day's  run  were 
turned  into  the  University  Computer  Center  which  has  a 
remote  terminal,  batch  desk  reader  and  interactive 
terminals  located  at  the  Engineering  Research  Center.  Once 
a  tape  was  in  the  computer  center,  a  simple  and  quick 
program  was  run  to  verify  that  the  data  on  the  tape  were 
correct  in  format  and  that  the  tape  was  readable  by  the 
computer.  The  data  on  the  tape  were  then  transferred  to  a 
"master"  tape  contained  internally  at  the  computer  center. 
All  actual  data  analysis  was  done  from  the  master  tape. 

All  data  channels  were  sampled  at  a  rate  of  about  one  scan 
per  2.25  seconds  resulting  in  about  800  discreet  points  per 
data  channel  for  each  30  minute  test  run.  The  advantage  of 
computer  analysis  in  this  case  is  obvious  in  the  quantity 
of  data  processed  during  analysis.  The  short  time  interval 
between  data  samples  thus  permitted  close  tracking  of  rapid 
temperature  changes  in  various  parts  of  the  system. 
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Calculations  required  knowledge  of  mass  flow  rates, 
temperatures  and  temperature  differences  across  the 
furnace.  Thermocouple  and  thermopile  outputs  were  recorded 
in  order  to  determine  the  heat  impact  from  the  furnace. 
Thermocouples  were  recorded  directly  in  absolute 
temperature  relative  to  the  ice  point  of  water  and  the 
thermopiles  were  recorded  as  a  voltage. 

The  thermopile  voltage,  or  generated  EMF  is  a  function  of 
the  difference  in  temperature  between  the  two  ends  of  the 
thermopile  and  the  number  of  junctions  making  up  the 
thermopile.  Because  thermocouples  are  not  linear  in  output 
with  temperature,  the  procedure  for  calculating  the  AT 
being  measured  requires  a  knowledge  of  the  absolute 
temperature  at  one  end  of  the  thermopile.  This  temperature 
becomes  the  thermopile  reference  temperature  for  a  series 
of  calculations  to  determine  the  at  at  the  thermopile. 

National  Bureau  of  Standardsl  equation  describing  the 
characteristics  of  T-type  thermocouples  were  used  in  these 
calculations.  The  reference  temperature  is  first  converted 
to  an  equivalent  EMF  relative  to  32°F.  Assuming  that  the 
reference  thermocouple  is  reading  the  low  temperature  side 
of  the  thermopile,  this  EMF  is  then  added  to  the  thermopile 
voltage  divided  by  the  number  of  thermocouple  junctions  in 
the  thermopile. 


EMF  =  EMF  T-PILE  VOLTAGE 

T-PILE  REFERENCE  +  N°  JUNCTIONS 


EMFT_piie  is  now  treated  as  though  it  were  itself  a 
thermocouple  voltage  and  using  the  NBS  equation  which 
relates  Temperature  to  EMF  the  high  side  thermopile 
temperature  is  calculated.  Then  At  is  this  high  side 
calculated  temperature  minus  the  original  reference 
temperature. 

The  Datametrics  Barocel  differential  pressure  sensor  was 
calibrated  at  the  Engineering  Research  Center 
Instrumentation  Shop  just  prior  to  the  furnace  testing. 


1  NBS  Monograph  125  (Thermocouple  Reference  Tables  Based 
on  the  IPTS-  68),  U.  S.  Department  of  Commerce,  National 
Bureau  of  Standards,  March,  1974. 
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The  output  voltage  of  this  unit  is  related  to  a 
differential  pressure  across  whatever  flow  device  is  being 
selected  by  both  a  sconxvalue  multi-port  switch  and  both 
the  switch  position  and  the  P  voltage  are  recorded. 

From  this  information,  mass  flow  rates  in  the  system  are 
calculated  in  order  to  finally  arrive  at  heat  delivery 
calculations.  In  the  computer  program,  these  calculations 
are  carried  out  for  every  2.25  second  increment  and  summed 
up  for  the  total  furnace  output. 

Furnace  thermal  efficiency  was  defined  as  the  heat 
delivered  to  the  room  circuit  divided  by  the  total  gas 
input  Btu's. 

HEAT  DELIVERED  TQ  ROOM.  BTU  X  10Q 
PERCENT  EFF  =  GAS  INPUT,  BTU 

The  jacket  loss  and  direct  losses  in  the  furnace  room  were 
not  considered  as  usable  heat  in  the  calculations. 

All  flow  rates  and  heating  values  were  converted  to 
standard  conditions  prior  to  calculating  efficiency.  The 
gas  heating  value  was  determined  from  gas  chromatograph 
analysis  compared  to  a  reference  standard.  During  the 
testing  period,  the  heating  value  of  the  gas  varied  from  a 
low  ot  1051  Btu/Ft3  to  a  high  of  1059  Btu/Ft3.  The 

actual  heating  values  were  inserted  in  the  computer  program 
for  analyzing  the  data  even  though  this  was  a  very  small 
change  at  less  than  1  percent. 

An  energy  balance  was  performed  for  the  furnace  using  the 
data  available  in  the  computer  program.  As  a  check  on  the 
results  obtained.  Gas  input  =  Heat  delivered  to  rooms  + 

Heat  up  the  stack  +  Jacket  losses. 

The  heat  delivered  to  the  rooms  is  made  up  of  two 
components.  The  difference  in  air  flow  measured  by  the 
return  air  monitor  and  the  delivery  air  monitor  is  the 
leakage  component  into  the  system  at  the  blower  cabinet. 
This  value  was  typically  around  180  Ft3/min  and  this 
quantity  of  flow  is  raised  from  furnace  room  temperature  to 
room  air  delivery  temperature.  The  significantly  larger 
component  is  that  flow  delivered  directly  through  the  room 
circuit  and  this  flow  is  raised  from  room  return 
temperature  to  room  delivery  temperature  for  heat  flow 
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calculations.  Jacket  losses  were  taken  to  be  1  percent  ot 
the  value  of  all  other  heat  quantities  combined. 


Closure  of  the  energy  balance  varied  with  duty  cycle.  The 
biggest  error  in  closure  was  for  the  10  percent  duty  cycles 
and  was  typically  somewhere  between  12  to  38  percent  otf. 
for  the  longer  duty  cycles  of  75  and  100  percent,  closure 
was  usually  within  2  or  3  percent,  but  occasionally  as  high 
as  5  percent  off. 

The  output  from  the  computer  program  was  a  complete  table 
of  temperatures,  flows,  calculated  temperature  differences 
and  indication  of  time  into  cycle  and  thermostat  and  blower 
status  indication.  Additionally,  the  program  would  output 
selected  temperature  plots  displaying  temperature  changes 
through  the  entire  30  minute  duty  cycle.  Average  outdoor 
temperature  was  calculated  and  printed  as  were  results  of 
the  energy  balance  along  with  the  individual  values  used  in 
the  balance  calculation. 


Results 

Appendix  D  contains  a  complete  description  of  the  furnace 
simulation  program  used  to  determine  the  seasonal  effects 
of  outside  combustion  air  on  furnace  performance.  Table  1 
is  a  compilation  of  the  key  results  from  the  simulation 
program  for  Boston,  MA;  Great  Falls,  MT;  and  Albuquerque, 
NM.  Two  types  of  thermostats  were  used  in  the  simulation, 
one  which  kept  the  house  at  a  constant  70of  24  hours  a 
day,  and  one  which  included  a  setback  to  60o  from  11  PM 
to  7  AM  each  night.  In  every  case,  the  outside  combustion 
air  reduced  the  seasonal  efficiency  calculated  by  the 
simulation  program  (modeled  seasonal  efficiency).  Fuel 
Ufifi.,  however,  was  reduced  by  outside  combustion  air 
because  of  the  reduction  in  templered  air  use  during 
furnace  on  cycles  and  the  reduction  in  tempered  air  loss  up 
the  stack  during  furnace  off  cycles. 
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TABLE  l 


ANNUAL  FUEL  USAGE  AND  MODELED  SEASONAL  EFFICIENCY 


Inside  Combustion  Air 

Outside  Combustion  Air 

Locat Ion 

Thennoacat 

Type 

Gas 

Input  , 

B  tux  10 

Room  Air 
Delivery 
BtuxlO*6 

Modeled 

Seasonal 

Efficiency 

Z 

Gas 

Input  . 
BtuxlO 

Room  Air 
Delivery 
BtuxlO 

Modeled 

Seasonal 

Efficiency 

Z 

Albuquerque, 

NM 

No  Night 
Setback 

187.5 

131.4 

70.1 

180.3 

124.0 

68.8 

Night 

Setback 

161.6 

114.0 

70.5 

155.7 

107.8 

69.2 

Bos  Con,  MA 

No  Night 
Setback 

192.6 

134.3 

69.7 

183.2 

125.0 

68.2 

Night 

Set  back 

169.6 

119.2 

70.3 

161.4 

111.1 

68.8 

ireat  Falls, 

NT 

No  Night 
Setback 

163.8 

113.7 

69.4 

153.9 

103.1 

67.0 

Night 

Setback 

146.6 

103.0 

70.3 

137.9 

93.5 

67.8 

Comparing  just  the  results  for  no  setbacks  on  the 
thermostat,  the  furnace  located  in  Albuquerque  ran  at  an 
efficiency  of  70.1  percent  and  used  187.5  million  Btu  of 
gas  to  heat  the  modeled  house  with  inside  combustion  air. 
With  outside  combustion  air,  the  efficiency  of  the  furnace 
dropped  slightly  to  68.8  percent  but  the  gas  input  needed 
to  heat  the  house  also  droplped  to  180.3  million  Btu. 
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Using  Boston  weather  tapes,  the  same  trend  was  indicated 
with  the  furnace  having  a  69.7  percent  efficiency  compared 
to  68.2  percent  and  a  gas  input  of  192.6  million  Btu 
compared  to  183.2  million  Btu  for  inside  compared  to 
outside  combustion  air  respectively.  The  same  results  for 
Great  Falls,  Mt,  were  69.4  percent  compared  to  67.0  percent 
for  efficiency  and  163.8  million  Btu  compared  to  153.9 
million  Btu  for  inside  compared  to  outside  combustion  air. 


The  reduction  in  gas  consumption  as  the  weather  conditions 
become  more  severe  is  the  result  of  sizing  the  house  load 
to  the  furnace  model  for  each  location.  Table  2  shows  the 
design  parameters  used  at  each  location  as  well  as  the 
actual  minimum  low  temperature  obtained  on  the  weather 
tapes. 

For  Albuquerque,  the  minimum  design  temperature  was 
exceeded  at  least  once  while  for  Boston,  it  was  equalled  at 
least  once.  For  Great  Falls,  MT,  however,  the  minimum 
design  temperature  was  never  exceeded  or  equalled.  This, 
combined  with  sizing  the  house  load  at  each  location  to  the 
furnace,  accounts  for  what  might  first  appear  to  be  an 
incorrect  trend  for  gas  consumption  in  the  simulation 
results. 

Table  3  presents  information  on  the  savings  achieved  for 
each  of  four  different  comparisons.  Concentrating  only  on 
the  comparison  of  indoor  to  outdoor  combustion  air  with  no 
setback  on  the  thermostat,  a  3.8  percent  savings  is 
achieved  for  Albuquerque  while  for  Boston  the  savings  are 
4.88  percent  and  for  Great  Falls,  MT,  they  are  6.04 
percent.  This  same  ranking  applied  to  the  comparison  of 
inside  and  outside  combustion  air  with  a  setback  thermostat 
although  the  savings  are  slightly  reduced. 

As  a  matter  of  interest,  the  savings  comparisons  for  the 
setback  and  nonsetback  thermostat  were  compared  and  can  be 
found  as  Columns  C  and  D  in  Table  3.  These  savings  show  a 
difference  from  the  comparisons  to  outdoor  and  indoor  air. 
That  is,  the  largest  percent  savings  now  occur  in 
Albuquerque  and  the  smallest  percent  savings  take  place  in 
Great  Falls,  MT. 

Table  4  snows  the  actual  annual  savings  calculated  from 
Table  3  using  the  October  1981  consumer  gas  price  in  each 
of  the  three  locations.  Annual  savings  attributable  to 
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inside  versus  outside  air  were  in  the  $20-50  range, 
suggests  that  payback  periods  for  retrofit  would  be 
over  ten  years. 


This 

well 


TABLE  2 


LOAD  DESIGN  PARAMETERS 


Design 

Temperature 

Actual 

T0Mtn 

Load 

UA 

BtuH/  F 

14.0 

8.9 

1050 

6.0 

6.0 

920 

-20.0 

-17.0 

650 

TABLE  3 


COMPUTED  SAVINGS  FOR  THE  VARIOUS  TEST  CASES 
SAVINGS  % 


Location 

A 

» 

C 

D 

Albuquerque 

1.84 

3.65 

13.81 

13.64 

Boston 

4.88 

4 .83 

n.94 

n.9o 

Great  fall* 

6.04 

5.93 

10.50 

10.40 

». 


Inside  vs.  Outside  Combustion  Air  -  Ho  Setback  Thermostat 
Ins ids  vs .  Outside  Combustion  Air  *  With  Setback  Thermostat 
Setback  vs.  No  Setback  Thermostat  -  Inaid®  Air 
S« check  vs.  No  Setback  Thermostat  -  (fetside  Bir 


(Inside -Outside) 

(  Inside  ) 

( Ins ide^Out side ) 

(  Inside  ) 

(No  Sstbaefc-Setback) 
(  No  Setback  ) 

(Ho  Setback-Setback) 
(  No  Setback  ) 
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TABLE  4 


DOLLAR  SAVINGS  FOR  THE  VARIOUS 
TEST  CASES,  $  PER  YEAR 


Location 

A 

B 

C 

D 

E 

Albuquerque 

$28.80 

$23.60 

$103.60 

$  98.40 

$127.20 

Boston 

54.05 

47.15 

132.25 

125 . 35 

179.40 

Great  Falls 

48.51 

42.63 

84.28 

78.40 

126.91 

A.  Inside  vs.  Outside  Combustion  Air  -  No  Setback  Thermostat 

(Inside-Outside) 

(  Inside  ) 

B.  Inside  vs.  Outside  Combustion  Air  -  With  Setback  Thermostat 

(Inside-Outside) 

(  Inside  ) 

C.  Setback  vs.  No  Setback  Thermostat  -  Inside  Air 

(No  Setback-Setback) 

(  No  Setback  ) 

D.  Setback  vs.  No  Setback  Thermostat  -  Outside  Air 

(No  Setback-Setback) 

(  No  Setback  ) 

E.  No  Setback,  Inside  Air  vs.  Setback,  Outside  Air 

(No  Setback,  Inside  -  Setback,  Outside) 

(  No  Setback,  Inside  ) 
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Table  5  is  a  listing  of  all  furnace  cycle  thermal 
efficiencies  measured  during  the  testing  program.  The 
values  in  the  table  are  listed  under  column  headings  of  the 
appropriate  outdoor  air  temperature  under  which  the  test 
runs  fall.  Listed  with  each  value  of  efficiency  is  the 
actual  average  outdoor  air  temperature  during  the 
individual  test  runs.  These  values  vary  somewhat  from  the 
actual  column  heading  but  are  grouped  under  the  nearest 
temperature  heading  for  placement  in  a  table  format.  Gaps 
in  data  shown  for  the  higher  temperatures  are  typically  the 
result  of  test  runs  presumed  to  be  that  temperature  at  the 
time  the  test  was  run,  but  on  analysis  of  the  data,  to  be 
closer  to  an  adjacent  column  heading. 

Testing  began  26  January  1981  and  continued  through  early 
April  1981  with  data  acquisition  proceeding  as  described  in 
Appendix  A.  Because  ot  the  need  to  obtain  cold  outdoor  air 
temperatures  as  they  became  available  at  this  time  of  year, 
testing  began  immediately  with  the  outdoor  vents  open 
rather  than  testing  of  the  so  called  "baseline"  or  indoor 
vent  condition. 

There  was  an  unusually  mild  winter  through  the  Rocky 
Mountain  Region  and  the  extreme  cold  temperatures  (zero 
degrees  F  and  below)  were  not  obtainable  until  a  single 
cold  period  on  the  10th  and  11th  of  February.  Testing 
proceeded  24  hours  a  day  in  order  to  ootain  all  data 
possible  during  this  cold  snap.  These  were  the  only  two 
days  of  near  zero  or  sub-zero  temperatures  for  the  entire 
winter.  Obviously,  data  for  test  runs  which  were  not 
obtained  during  this  two-day  period  were  not  available. 

Figures  1  through  6  are  plots  of  the  data  from  Table  5. 

Each  duty  cycle  from  10  percent  to  100  percent  is 
represented  as  a  separate  plot  with  thermal  efficiency 
plotted  against  outdoor  air  temperature  for  all  wind 
conditions. 

Wind  appears  to  have  little  or  no  effect  on  furnace 
efficiencies.  Points  representing  moderate  and  high  wind, 
as  shown  in  Table  6,  for  given  temperatures,  fall  within 
the  same  range  of  efficiencies  for  no  wind.  Later,  in 
examining  plots  of  actual  temperatures  in  the  system,  there 
is  noted  an  effect  on  stack  temperatures.  The  effect  of 
wind  is  apparently  confined  primarily  to  dilution  air  at 
the  draft  hood  and  does  not  alter  the  actual  combustion  or 
efficiency  or  the  heat  exchanger.  Table  2  shows  the 
measured  effect  of  wind  on  stack  flow. 
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Cold  outside  combustion  air  has  the  most  dramatic  effect  on 
furnace  efficiencies  for  the  short  duty  cycles.  Figure  1, 
the  plot  of  efficiency  versus  outdoor  temperature  for  the 
10  percent  duty  cycle,  shows  a  range  of  efficiencies  from 
about  42  percent  for  a  -5°F  outdoor  temperature  to  a  high 
of  about  58  percent  for  60°F  outside  temperature.  There 
is  considerable  scatter  in  the  data  for  the  10  percent  test 
as  also  indicated  by  the  largest  error  in  the  energy 
closure  equations  for  this  duty  cycle.  The  general  shape 
and  slope  of  the  curve  fit  to  these  points  is  carried 
forward  from  the  plots  of  the  longer  duty  cycles  where 
scatter  is  considerably  less.  The  10  percent  duty  cycle 
shows  a  total  change  in  efficiency  of  16  percent  from  - 
5of  to  60of. 

As  the  "burner  on"  duty  cycle  increases,  the  low  and  high 
values  of  efficiency  for  the  same  temperatures,  -5of  and 
60°F,  gradually  improve  until  the  best  efficiencies  are 
obtained  at  a  100  percent  duty  cycle  ,  Figure  6.  At  -50f 
outside  air  temperature,  the  furnace  efficiency  for  the 
100  percent  duty  cycle  is  approximately  69  percent  and  this 
improves  by  only  3  percent  to  a  value  of  72  percent  for 
60°F  outside  air. 
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HIGH  WIND  MODERATE  WIND  NO  WIND 


TABLE  5 


FURNACE  THERMAL  EFFICIENCIES 


Combustion  &  Dilution  Air 


Outside  Air,  Temperature  of 


Percenrv 

Duty  Cyclek 

1 

-10 

0 

10 

25 

A0 

55 

10 

__ 

45.5 

48.4 

60.9 

57.4 

57.7 

(-0.1) 

(12.1) 

(21.9) 

(47.2) 

(53.6) 

20 

— 

56.4 

55.5 

60.8 

64.5 

60.7 

(-1.3) 

(11.7) 

(24.0) 

(45.3) 

(58.7) 

§  30 

61.2 

62.1 

61.2 

67.3 

64.0 

66.0 

M 

5 

(-4.2) 

(-3.0) 

(11.1) 

(22.6) 

(36.7) 

(57.3) 

o  50 

65.2 

66.7 

67.1 

69.5 

68.9 

69.3 

z 

(-4.6) 

(-2.9) 

(11.7) 

(26.7) 

(45.7) 

(60.1) 

75 

67.4 

68.4 

69.6 

72.4 

72.6 

— 

(-6.8) 

(-3.5) 

(15.1) 

(26.7) 

(41.5) 

100 

66.6 

70.0 

67.0 

69.2 

— 

72.2 

(-6.9) 

(0.3) 

(18.7) 

(27.2) 

(55.8) 

10 

— 

_ 

— 

_ 

59.2 

52.8 

(44.0) 

(56.7) 

a  20 

— 

— 

— 

— 

64.4 
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FIGURE  1.  Plot  of  Cyclic  Thermal  Efficiency  Versus  Outdoor  Temperature  for 
10  Percent  Duty  Cycle,  Using  Outdoor  Combustion  and  Dilution 
Air.  Difference  in  Efficiencies  is  16  Percent,  From  -5°  to  60°. 


OUTDOOR  TEMPERATURE,  °F 


FIGURE  2.  Plot  of  Cyclic  Thermal  Efficiency  Versus  Outdoor  Temperature  for 
20  Percent  Duty  Cycle,  Using  Outdoor  Combustion  and  Dilution 
Air.  Difference  in  Efficiencies  is  9  Percent,  From  -5°  to  60°. 
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FIGURE  3.  Plot  of  Cyclic  Thermal  Efficiency  Versus  Outdoor  Temperature  for 
30  Percent  Duty  Cycle,  Using  Outdoor  Combustion  and  Dilution 
Air.  Difference  in  Efficiencies  is  7  Percent,  From  -5°  to  60°. 


FIGURE  4.  Plot  of  Cyclic  Thermal  Efficiency  Versus  Outdoor  Temperature  for 
50  Percent  Duty  Cycle,  Using  Outdoor  Combustion  and  Dilution 
Air.  Difference  in  Efficiencies  is  6  Percent,  From  -5°  to  60°. 


OUTDOOR  TEMPERATURE,  °F 

FIGURE  5.  Plot  of  Cyclic  Thermal  Efficiency  Versus  Outdoor  Temperature  for 
75  Percent  Duty  Cycle,  Using  Outdoor  Combustion  and  Dilution 
Air.  Difference  in  Efficiencies  is  6  Percent,  From  -5°  to  60°. 
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FIGURE  6.  Plot  of  Cyclic  Thermal  Efficiency  Versus  Outdoor  Temperature  for 
100  Percent  Duty  Cycle,  Using  Outdoor  Combustion  and  Dilution 
Air.  Difference  in  Efficiencies  is  3  Percent,  From  -5°  to  60°. 
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TABLE  6 


WIND  INDUCED  DRAFT  EFFECT 


Burner  Off 
Burner  On 


Vent  Stack  Flow  Rates 


25  MPH  Wind  40-50  MPH  Wind 


No  Wind 
20  ft^/min. 
65  ft^/min. 


55  ft^/min. 

90  ft'Vmin. 


70  ft'Vmin. 
100  ftVmin. 


The  change  in  efficiency  over  the  -5of  to  60of 
temperature  range  is  highest  for  the  10  percent  duty  cycle 
and  gradually  decreases  through  longer  duty  cycles;  16 
percent  change  at  10  percent  duty  cyle,  9  percent  change  at 
20  percent  duty  cycle,  7  percent  change  at  30  percent  duty 
cycle,  6  percent  change  for  both  50  percent  and  75  percent 
duty  cycles,  and  finally,  only  a  3  percent  change  at  100 
percent  duty  cycle.  At  the  same  time,  the  actual 
efficiencies  covered  by  the  ranges  gradually  increase  with 
longer  duty  cycles;  42  to  58  percent  for  10  percent  duty 
cycle,  55  to  69  percent  for  20  percent  duty  cycle,  61  to 
68  percent  for  30  percent  duty  cycle,  65  to  71  percent  for 
50  percent  duty  cycle,  66  to  72  percent  for  75  percent  duty 
cycle,  and  69  to  72  percent  for  100  percent  duty  cycle. 

Figures  7  through  12  are  plots  of  the  same  test  runs  as 
Figures  1  through  6  but  instead  of  outdoor  air  temperature, 
furnace  efficiency  is  plotted  versus  furnace  room 
temperature.  This  furnace  room  temperature  is  the  average 
temperature  during  the  run,  but  even  at  the  100  percent 
duty  cycles,  it  did  not  vary  appreciably.  Also,  the 
average  furnace  room  temperature  at  a  given  time  was  the 
average  of  5  sensors  located  around  the  furnace  room. 

These  figures  do  not  display  any  new  information  even 
though  the  real  variable  affecting  furnace  performance 
would  be  the  actual  air  temperature  in  the  furnace  room. 

The  plots  are  presented,  therefore,  as  a  matter  of 
interest . 
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FIGURE  7.  Plot  of  Cyclic  Thermal  Efficiency  Versus  Furnace  Room  Tempera¬ 
ture  for  10  Percent  Duty  Cycle,  Using  Indoor  and  Outdoor 
Combustion  and  Dilution  Air.  Difference  in  Efficiencies  is  19 
Percent,  From  20°  to  80°. 


FIGURE  8.  Plot  of  Cyclic  Thermal  Efficiency  Versus  Furnace  Room  Tempera¬ 
ture  for  20  Percent  Duty  Cycle,  Using  Indoor  and  Outdoor 
Combustion  and  Dilution  Air.  Difference  in  Efficiencies  is  12 
Percent,  From  20°  to  80°. 
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Plot  of  Cyclic  Thermal  Efficiency  Versus  Furnace  Room  Tempera¬ 
ture  for  10  Percent  Duty  Cycle,  Using  Indoor  and  Outdoor 
Combustion  and  Dilution  Air.  Difference  in  Efficiencies  is  10 
Percent,  From  20°  to  80°. 
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FIGURE  10.  Plot  of  Cyclic  Thermal  Efficiency  Versus  Furnace  Room  Tempera¬ 
ture  for  50  Percent  Duty  Cycle,  Using  Indoor  and  Outdoor 
Combustion  and  Dilution  Air.  Difference  in  Efficiencies  is  6 
Percent,  From  20°  to  80°. 
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FIGURE  11*.  Plot  of  Cyclic  Thermal  Efficiency  Versus  Furnace  Room  Tempera¬ 
ture  for  75  Percent  Duty  Cycle,  Using  Indoor  and  Outdoor 
Combustion  and  Dilution  Air.  Difference  in  Efficiencies  is  7 
Percent,  From  20°  to  80°. 
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FIGURE  12.  Plot  of  Cyclic  Thermal  Efficiency  Versus  Furnace  Room  Tempera¬ 
ture  for  100  Percent  Duty  Cycle,  Using  Indoor  and  Outdoor 
Combustion  and  Dilution  Air.  Difference  in  Efficiencies  is  3 
Percent,  From  20°  to  80°. 
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Temperature  profiles  for  the  30-minute  test  cycle  for 
selected  test  runs  are  shown  as  Figures  13  through  22. 

These  plots  are  shown  for  various  test  duty  cycles  and  for 
various  outdoor  temperatures  and  wind  conditions. 

Figures  13  and  14  are  temperature  profile  plots  for  a  10 
percent  duty  cycle  with  outside  air.  Figure  13  is  for 
still  air  conditions  while  Figure  14  is  for  a  high  wind 
condition.  With  the  exception  of  only  one  recorded 
temperature,  the  plots  are  essentially  similar. 

Temperatures  in  the  air  delivery  duct,  the  furnace  plenum 
and  the  furnace  room  are  essentially  the  same  from  plot  to 
plot.  Only  the  vent  stack  temperature  shows  any 
appreciable  effect  from  the  cold  outside  combustion  air 
combined  with  wind  effect.  This  temperature  has  dropped 
nearly  65°F  from  a  value  of  265of  with  no  wind  to 
200OF  with  a  high  wind  blowing. 

Figures  15  and  16  are  plots  of  a  10  percent  duty  cycle  with 
the  combustion  air  source  as  inside  tempered  air.  Figure 
15  is  a  still  air  condition,  and  Figure  16  is  a  high  wind 
condition.  Almost  the  same  situation  takes  place  as  with 
the  outside  combustion  air  except  the  wind  condition  stack 
temperature  is  not  quite  as  low.  In  this  case,  the  stack 
temperature  with  no  wind  is  about  260°F  which  then  drops 
about  50°F  to  a  value  ot  210<>F  with  a  high  wind.  The 

furnace  room  temperatures  are  substantially  different  for 
the  outside  compared  to  the  inside  combustion  air  which 
accounts  for  the  change  in  stack  temperature  lows  in  the 
high  wind  case. 

Figures  17  and  18  demonstrate  the  same  effect  for  a 
slightly  longer  duty  cycle  of  30  percent.  In  this  case, 
the  wind  is  a  moderate  one  and  the  effect  registers  only  on 
stack  temperature.  The  other  temperatures  plotted  in  the 
figures  are  virtually  identical  in  the  two  tests  and  the 
computed  furnace  efficiencies  are  within  2  percent  at  62.1 
percent  and  59.8  percent,  respectively. 


FURNACE  ROON 


FIGURE  13.  Plot  of  Temperatures  Versus  Time  into  Cycle  for  10  Percent 
Duty  Cycle  With  No  Wind,  Using  Outside  Combustion  and 
Dilution  Air.  Furnace  Efficiency  is  60.9  Percent. 


FIGURE  14.  Plot  of  Temperatures  Versus  Time  into  Cycle  for  10  Percent 
Duty  Cycle  With  High  Wind,  Using  Outside  Combustion  and 
Dilution  Air.  Furnace  Efficiency  is  41.8  Percent. 


Duty  Cycle  With  No  Wind,  Using  Outside  Combustion 
Dilution  Air.  Furnace  Efficiency  is  61.2  Percent. 


figure  18.  Plot  of  Temperatures  Versus  Time  into  Cycle  for  30  Percent 

Duty  Cycle  With  Moderate  Wind,  Using  Outside  Combustion  and 
Dilution  Air.  Furnace  Efficiency  is  59.8  Percent. 
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FIGURE  19.  Plot  of  Temperatures  Versus  Time  into  Cycle  for  75  Percent 
Duty  Cycle  With  No  Wind,  Using  Outside  Combustion  and 
Dilution  Air.  Furnace  Efficiency  is  67.4  Percent. 
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FIGURE  20.  Plot  of  Temperatures  Versus  Time  into  Cycle  for  75  Percent 
Duty  Cycle  With  Moderate  Wind,  Using  Outside  Combustion  an 
Dilution  Air.  Furnace  Efficiency  is  68.3  Percent. 
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FIGURE  21.  Plot  of  Temperatures  Versus  Time  into  Cycle  for  75  Percent 
Duty  Cycle  With  High  Wind,  Using  Outside  Combustion  and 
Dilution  Air.  Furnace  Efficiency  is  66.7  Percent. 
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FIGURE  22.  Plot  of  Temperatures  Versus  Time  into  Cycle  for  75  Percent 
Duty  Cycle  With  No  Wind,  Using  Inside  Combustion  and 
Dilution  Air.  Furnace  Efficiency  is  72.6  Percent. 


Figures  19,  20  and  21  are  75  percent  duty  cycle  runs  with 
outside  combustion  air  at  Oof  or  below.  The  furnace 
efficiencies  for  these  test  runs  were  67.9  percent,  68.3 
percent  and  66.7  percent,  respectively  for  still  air, 
moderate  wind  and  high  wind.  Again,  no  significant  effect 
on  furnace  efficiency  is  indicated  with  only  the  vent  stack, 
temperature  showing  any  change  from  plot  to  plot.  Figure 
22  is  a  plot  of  inside  combustion  air  and  no  wind  for  a  75 
percent  duty  cycle  for  comparison  to  the  previous  three 
figures.  The  furnace  efficiency  for  this  test  run  was 
72.6  percent  and  the  furnace  room  temperature  was  about 
70OF. 


CONCLUSIONS 

Furnace  cyclic,  thermal  efficiency  decreased  with  colder 
outdoor  temperatures  when  outdoor  combustion  and  dilution 
air  was'  brought  into  the  furnace  room. 

The  decrease  in  turnace  efficiency  with  decreasing  outdoor 
temperature  was  most  pronounced  for  the  10  percent  duty 
cycle  wnere  efficiency  changed  by  16  percent  within  the 
tested  temperature  range.  There  was  significantly  less 
drop  in  efficiency  with  decreasing  outdoor  temperature  as 
the  duty  cycle  increased  until  at  100  percent  duty  cycle, 
the  change  in  efficiency  was  only  3  percent  over  the 
temperature  range  tested. 

When  outdoor  combustion  air  was  being  used,  temperatures 
within  the  furnace  room  dropped  below  freezing  even  during 
the  100  percent  duty  cycle.  This  could  present  problems  in 
locating  a  gas-fired  water  heater  in  the  furnace  room, 
sharing  a  common  vent,  because  of  possible  freezing.  Also 
losses  from  the  water  tank  to  the  cold  surrounding  air  may 
be  substantially  increased  unless  additional  insulation  is 
applied  to  the  tank. 

Wind-induced  draft  had  no  effect  on  furnace  cyclic,  thermal 
efficiency  even  when  the  temperature  in  the  furnace  room 
was  quite  low. 

Computer  simulations  run  for  a  furnace  model  and  house  at 
three  locations  in  the  U.S.  showed  that  outside  combustion 
air  resulted  in  an  annual  fuel  savings,  ranging  from  nearly 
4  percent  to  a  little  better  than  6  percent,  when  compared 
to  a  furnace  and  nouse  using  inside  combustion  air. 


The  savings  achieved  by  using  outside  combustion  air  were 
largest  for  the  coldest  climate  tested  and  smallest  for  the 
more  moderate  climate. 


RECOMMENDATIONS 

Retrofit  of  residential  gas  furnaces  should  be  considered 
in  all  climates  based  upon  a  cost-effectiveness  study  of 
each  furnace  circumstance. 

Cost  of  retrofit  should  be  calculated  based  upon  compliance 
with  the  local  building/mechanical  codes  and  NFPA  54,  the 
National  Fuel  Gas  Code. 

Consider  automatic  vent  dampers  as  an  alternative  to 
outside  combustion  air  and  compare  these  two  candidate 
retrofits.  Refer  to  Report  FESA-TS-2072 . 


41 


APPENDIX  A 

THE  TEST  FACILITY  AND  INSTRUMENTATION 
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The  test  facility  consisted  of  a  "furnace  room"  6  feet  wide 
and  12  feet  long  by  8  feet  high  and  an  instrumentation  room 
the  same  size  located  adjacent  to  the  furnace  room.  The 
test  rooms  were  located  inside  the  hydro-machinery 
laboratory  space  at  the  Engineering  Research  Center  at  the 
Foothills  Campus  of  Colorado  State  University.  The 
laboratory  has  a  large  enclosed  floor  of  90  feet  by  175 
feet.  This  floor  space  is  heated  in  the  winter  by  radiant 
style  heaters  at  the  root  level  with  vent  fans  venting  in 
the  roof. 

This  facility  was  still  intact  from  the  Task  Order  No.  7 
study  of  vent  energy-saving  devices  and  only  minor  changes 
and  additions  were  made.  As  originally  constructed,  the 
furnace  room  was  insulated  from  the  attached  instrumenta¬ 
tion  room  and  trom  the  surrounding  hydro-machinery  lab. 

The  concrete  exterior  wall  of  the  furnace  room  was  not 
insulated  nor  was  the  ceiling  space. 

No  special  attempt  was  made  to  seal  the  walls  or  doorway 
against  air  intrusion.  Duct  penetrations,  holes  for 
instrumentation  leads,  etc.  into  the  furnace  room  were 
sealed  to  form  a  reasonably  airtight  penetration.  This  is 
assumed  to  be  typical  of  how  a  retrofit  for  outdoor 
combustion  air  might  be  accomplished. 

The  "room  air"  for  the  test  furnace  was  supplied  from  the 
laboratory  space.  Heated  air  was  delivered  to  the  lab 
space  at  8  feet  above  the  floor  and  directed  in  an  opposite 
direction  from  the  cold  air  intake.  The  air  temperature 
near  the  floor  of  the  laboratory  varied  between  65°F  to 
70°F  which  is  representative  of  return  air  temperature  in 
a  heated  building. 

A  regulated  natural  gas  supply  line  already  existed  within 
the  building  and  was  easily  extended  to  the  test  rooms  for 
operation  of  the  furnace.  The  combustion  vent  was  run 
through  an  office  balcony  located  immediately  above  the 
test  rooms  and  vented  into  a  space  near  the  "roof"  inside 
the  building.  A  vent  fan,  operating  whenever  the  furnace 
was  operating,  then  pulled  the  combustion  gases  outside  the 
lab  building.  This  arrangement  enabled  tests  to  be 
independent  of  random  wind  effects  on  the  vent.  To 
simulate  a  wind  induced  draft  through  the  vent,  a  two-speed 
vent  fan  was  used.  This  allowed  the  effects  of  a  wind 
induced  draft  to  be  duplicated  consistently  from  run  to  run 
without  being  dependent  on  the  actual  wind.  The  vent  fan 
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was  installed  so  that  the  flue  gas  flow  would  not  be 
obstructed  when  the  fan  was  not  being  used. 

Combustion  and  dilution  air  for  furnace  operation  was 
provided  through  one  of  two  alternate  sets  of  combustion 
air  vents.  For  the  baseline  case,  where  combustion  and 
dilution  air  are  defined  as  tempered  indoor  air,  two  vents 
were  placed  in  a  wall  common  to  the  hydro-machinery  lab 
space.  These  vents  had  a  gross  area  of  144  in2  each  with 
one  located  with  center  11  inches  above  the  floor  and  the 
second  with  center  located  16  inches  below  the  ceiling. 

The  gross  area  was  reduced  to  an  effective  free  area  of  100 
in2  each  by  installation  of  a  metal  grill  over  the 
openings. 

The  outside  air  source  for  combustion  air  was  provided  by 
drilling  four  6  inch  diameter  holes  through  the  concrete 
wall,  two  located  with  centers  17  inches  above  the  floor 
and  two  with  centers  located  12  inches  below  the  ceiling. 
This  provided  a  gross  area  of  56.5  in2  for  both  the  low 
and  high  vent.  This  was  reduced  to  an  effective  free  area 
of  about  35  in2  for  each  of  the  high  and  low  vents  after 
a  grill  and  screen  cover  was  installed. 

Each  of  the  combustion  air  vents  are  fitted  with  a  tight 
removable  cover  so  that  air  vents  were  blocked 
appropriately  for  testing  either  for  outdoor  air  or  for 
indoor  air.  The  outdoor  air  vents  had  direct  access  to  the 
outside  air  along  the  main  west  wall  of  the  hydro-machinery 
lab  building. 

The  furnace  used  for  the  tests  was  a  Lennox  Model  G1203-110 
with  a  nameplate  rating  at  110,000  BtuH  input.  The  furnace 
is  a  complete  standard  model  just  as  shipped  by  the 
manufacturer.  The  unit  included  a  four  speed  blower  and  a 
standing  gas  pilot  light  and  all  burner  controls,  fan 
switches  and  24  volt  control  wiring.  The  furnace  was 
purchased  from  a  local  distributor  and  delivered  in  the 
original  factory  carton.  Figure  Al  is  an  overall  view  of 
the  furnace  as  installed  with  connecting  duct  work. 

Furnace  installation  was  accomplished  by  a  local  heating 
contractor.  Blower  fan  settings  were  placed  on  the  "medium 
high"  speed,  which  in  many  cases  accommodated  air 
conditioning  systems  installed  in  conjunction  with  the 
furnace  unit.  Fan  thermostat  settings  were  adjusted  so 
that  the  fan  would  turn  off  at  around  80of.  This  setting 


was  made  to  assist  in  smoothing  out  blower  operation  time 
through  the  various  tests. 

Air  ducts  to  the  furnace  from  the  laboratory  space  were 
stubbed  through  the  instrumentation  room  and  through  the 
common  wall  ot  the  furnace  room.  All  connections  on  the 
duct  were  sealed  with  a  silicone  sealant  to  prevent  air 
leakage  from  or  air  infiltration  into  the  "room  air" 
system.  The  ducts  were  wrapped  with  3-1/2  inches  of  fiber 
glass  insulation  with  an  R  value  of  11.  Two  Cambridge 
Filter  Corporation  air  flow  monitors  Model  FMS-D-14  X  14, 
consisting  ot  pitot  rakes  and  static  pressure  probes  were 
installed,  one  in  the  cool  air  return  duct  and  one  in  the 
warm  air  delivery  duct.  The  duct  air  flow  monitors,  were 
calibrated  at  the  Solar  Energy  Applications  Laboratory 
using  an  ASME  flow  nozzle  as  the  calibration  standard. 

Each  monitor  consisted  of  total  and  static  pressure  probe 
arrays  which  were  separately  manifolded  to  an  outlet  port. 
The  differences  between  total  and  static  pressures  were 
measured  with  a  Barocel  electronic  manometer.  Average  air 
velocity  at  the  cross-section  of  the  monitor,  hence 
volumetric  air  flow  rate  and  flow  rate  coefficients  were 
determined  from  the  calibrations.  Figure  A2  shows  an  air 
flow  monitor  as  installed  in  the  ductwork.  Figure  A3  is  an 
overall  plan  view  of  the  test  rooms  and  furnace  layout. 


Instrumentation 

The  following  data  were  taken  for  the  analysis  of  each 
device  on  the  vent  stack: 

1.  Temperatures  throughout  the  system. 

2.  Temperature  differentials  across  the  furnace. 


45 


3.  Air  flow  rates  in  the  room  air  circuit 


4.  Vent  flow  rates  during  and  before  burner  operation. 

5.  Wet  and  dry  bulb  temperature  in  the  lab. 

6.  Gas  consumption  during  each  run. 

7.  Stack  gas  samples  for  CC>2,  O2  and  CO. 

8.  Natural  gas  samples  for  heating  value  analysis. 


Temperatures 

Copper-constantan  (Type  T)  thermocouples  were  used  to 
monitor  temperatures  throughout  the  system  during  testing. 
Figures  A4  and  A5  are  system  schematics  giving  relative 
positions  of  temperature  sensors.  Thermocouples, 
designated  by  TC  on  the  figures  were  used  to  measure 
absolute  temperatures. 

The  thermocouples  were  formed  from  type  T  thermocouple  wire 
manufactured  to  special  limits  of  purity.  For  the  vent 
positions  at  a  high  temperature,  teflon- insulated 
thermocouple  wire  was  used  and  a  ceramic  thermal  coating 
was  applied  over  the  exposed  thermocouple  to  protect  it 
from  possible  deterioration.  The  "special  limits  of  error" 
(3/8  percent)  wire  used  results  in  thermocouples  of  more 
uniform  and  predictable  characteristics.  Junctions  were 
welded  bead  and  measures  were  taken  to  insure  a  uniform  and 
small  size  of  the  weld  bead  to  reduce  response  time. 

Thermopiles,  which  are  made  up  of  individual  thermocouple 
pairs  wired  in  series,  were  used  to  measure  all  temperature 
differences.  The  thermopile  circuits  do  not  measure 
absolute  temperature  relative  to  some  reference  but  rather, 
provide  a  measurement  of  temperature  increase  or  decrease 
(delta)  across  the  circuit  points  relative  to  either  end  of 
the  circuit.  For  example,  a  sixteen  junction  thermopile 
was  used  to  measure  differences  in  air  temperature  across 
the  furnace,  designated  as  TP-11  in  Figure  A2.  This 
thermopile  detects  the  increase  in  temperature  across  the 
furnace  between  the  cool  air  in  and  the  warm  air  delivery. 
Each  of  the  sixteen  thermocouples  making  up  the  thermopile 
was  placed  at  the  centroid  of  an  equal  area  in  the  duct 
cross  section.  In  this  way,  the  thermopile  provides  an 
average  of  the  air  temperature  through  the  duct 
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cross-section.  The  thermopile  junctions  were  mounted  to  an 
open  wire  grid  in  the  duct  section  and  carefully  fitted 
along  the  grid  to  cover  a  minimum  of  pressure  drop  in  the 
air  flow.  The  thermocouple  absolute  temperature  junction 
was  also  placed  at  the  center  of  this  grid. 
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FIGURE  A5.  View  of  Instrumentation:  (1)  Doric  Digital  Acquisition  Systo 
(2)  Kennedy  Incremental  Magnetic  Tape  Recorder,  (3)  Furnace 
as  Seen  Through  the  Observation  Window. 
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Data  Acquisition  System 

A  digital  data  acquisition  system  with  magnetic  tape 
recorder  was  used  to  collect  the  data  over  the  30-minute 
duty  cycle  with  enough  data  points  to  make  significant 
energy  flow  calculations.  The  system  used  was  a  Doric 
Digitrend  Model  220  linked  with  a  Kennedy  Incremental 
Model  1600  digital  tape  recorder.  The  Doric  sampled  all 
temperatures  and  thermopiles  approximately  every 
1.8  seconds  in  a  continuous  scan  mode  and  the  data  was 
recorded,  via  the  Kennedy,  on  magnetic  tape.  The  Doric 
system  has  temperature  conservation  circuits  contained 
within  the  wire  connection  blocks  so  that  thermocouple 
junctions  can  input  directly  and  read  out  in  absolute 
values  relative  to  the  ice  point  of  water.  This  instrument 
was  factory  calibrated  and  then  checked  after  receipt  on 
site  immediately  prior  to  the  testing  program.  Figure  A6 
is  an  overall  view  of  the  instrumentation  roack  containing 
the  Doric  and  Kennedy  data  acquisition  system. 

Table  A1  is  a  listing  of  the  channels  sampled  on  the  Doric 
during  test  runs.  The  information  consists  of  either  a 
thermocouple  temperature  measurement  or  a  voltage  with  the 
channel  range  approximately  set  to  accommodate  the  maximum 
voltage  expected  to  achieve  the  best  resolution.  Other 
T  type  thermocouples  were  read  directly  in  temperature,  the 
thermopiles  produce  a  voltage  and  are  read  on  a  multivolt 
scale. 

Several  other  channels  of  data  are  also  a  voltage 
relationship,  including  the  air  flow  measurements 
throughout  the  system,  a  signal  to  indicate  which  airflow 
measurement  is  being  sampled,  the  gas  temperature  and  the 
blower  and  thermostat-on  indicators. 
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TABLE  AJL 


LIST  OF  DATA  CHANNELS  RECORDED  ON  THE  DORIC  DIGITAL  AQUISITION  SYSTEM 


Doric  Data 
Channel 

Label 

1 

TCI  00 

2 

TC101 

3 

TCI  02 

4 

TC200 

5 

TC201 

6 

TC202 

7 

TC300 

8 

TC301 

9 

TC302 

10 

TC303W 

11 

TC304N 

12 

TC305E 

13 

TC306W 

14 

TC400 

15 

TC401 

16 

TC402 

17 

TD100 

18 

TD101 

19 

TD200 

20 

TD201 

21 

TD202 

22 

DP  101 

23 

T250 

24 

SP1 

25 

SP2 

26 

SP3 

27 

SP4 

28 

THERM 

29 

BLOW 

30 

TCI  03 

Units 


it 

R 

n 

r 

n 

n 

R 

r 

H 

R 

fl 

R 

fl 

R 

R 

mV 

R 

H 

R 

R 

R 


Room  Air  Return 
Room  Air  Delivery 
Room  Air  Delivery 
Stack 
Stack 

Upper  Stack  Air 
Furnace  Room  Air 
Low  Interior  Vent 
High  Interior  Vent 
Furnace  Room  Air  (West) 
Furnace  Room  Air  (North) 
Furnace  Room  Air  (East) 
Furnace  Room  Air  (South) 
Outside  Air 
Low  Exterior  Vent 
High  Exterior  Vent 
Air  Circuit  Thermopile 
Air  Circuit  Thermopile 
Stack  Circuit  to  Room  Air 
Stack  Circuit  to  Room  Air 
Stack  Circuit 

Barocel  Differential  Pressure 
Gas  Temperature  (Thermistor) 
Scanivalve  Position  for 
Barocel 

(Binary  Coded) 

Thermostat  on  Indicator 
Blower  on  Indicator 
Furnace  Plenum 


FIGURE  a6.  View  of  Instrumentation  Rack.  From  the  Top:  (1)  Thermostat 
Timer,  (2)  Kennedy  Incremental  Magnetic  Tape  Recorder, 

(3)  Doric  Digital  Acquis* tion  System  (4)  Barocel  Differential 
Pressure  Meter. 


Air  Plow  Rates 


The  air  flow  rates  in  the  room  air  duct  system  during 
blower  operation  and  the  stack  flow  rates  are  required  to 
compute  the  energy  flux  from  the  furnace. 

Two  Cambridge  Filter  Corporation  duct  air  flow  monitors. 
Model  PMS-D  were  used  in  the  room  air  system:  one  in  the 
return  air  duct  and  one  in  the  air  delivery  duct.  These 
flow  monitors  were  located  to  have  a  long  straight  duct 
section  at  the  upstream  approach.  Although  care  was  taken 
to  seal  all  duct  connections  with  flexible  silicone 
sealant,  no  special  measures  were  taken  with  the  furnace 
blower  cabinet  and  the  two  flow  minitors  are  needed  to 
provide  the  air  leakage  rate  into  the  system. 

The  flow  monitors  consist  of  a  series  of  pitot  tubes,  known 
as  a  pitot  rake,  which  are  connected  into  a  common 
averaging  manifold.  These  individual  pitot  tubes  are 
arranged  into  a  pattern  that  samples  velocity  at  equal  area 
intervals  through  the  flow  monitor  cross-section.  Also 
contained  within  is  a  second  set  of  static  pressure  parts 
brought  out  to  a  second  averaging  manifold.  The  two 
manifolds  connect  to  ports  on  the  side  of  the  flow  monitor 
and  the  pressure  difference  detected  here  is  proportional 
to  air  flow.  The  pressure  differential  was  detected  by 
running  tubing  from  the  flow  monitor  to  a  "scanivalue" 
selection  switch  and  then  to  a  Datametrics  Corp.  Barocel, 
and  Electronic  Model  1173  Manometer  which  is  a  low-level 
differential  pressure  transducer.  The  Barocel  converts  the 
pressure  to  a  voltage  which  is  sensed  by  one  of  the  Doric 
data  channels.  The  scanivalue  selector  was  manually 
stepped  at  regular  intervals  during  test  runs  to  sample  all 
the  air  flow  rates  in  the  system. 

Stack  flow  rates  were  measured  with  a  pitot  static  tube  at 
the  centerline  of  the  stack  located  in  a  straight  section 
of  stack  passing  through  the  balcony  room  above  the  test 
room.  Pressure  differentials  on  the  pitot  static  tube  were 
also  read  on  the  Barocel  by  stepping  the  scanivalue. 
Attempts  to  use  a  hot  wire  anemometer  to  measure  stack  flow 
rates  proved  unsuccessful  despite  the  manufacturer's  claim 
that  the  unit  used  should  function  in  the  stack 
environment. 
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Wet  and  Dry  Bulb  Temperatures 


A  Taylor  Model  1328  sling  psychrometer  was  used  to  measure 
dry  bulb  and  wet  bulb  temperatures,  both  in  laboratory 
space  and  m  the  furnace  room.  Readings  were  taken  for 
each  test  run  and  manually  recorded  on  the  data  sheets. 


Blower  Power 

The  electric  energy  to  the  blower  was  measured  with  a  power 
watt  meter.  This  value  was  manually  entered  in  the  data 
sheets  for  each  run.  While  parasitic  power  should  be 
subtracted  from  the  heat  supplied  by  the  burner,  blower 
energy  consumed  was  not  large  enough  to  have  a  significant 
effect  on  the  comparisons  of  the  different  testing  condi¬ 
tions. 


Gas  Input 

Volume  of  gas  input  to  the  test  furnace  was  measured  using 
a  Sprague  Corp  Model  175  Test  Meter  Serial  3433436 
specially  factory-cal ibrated  for  testing  purposes.  The 
meter  was  located  within  the  furnace  room  space  close  to 
where  the  gas  temperature  was  monitored.  The  meter 
readings  were  taken  manually  and  recorded  on  data  sheets 
during  the  test  runs. 

Gas  line  pressure  was  monitored  using  a  Dwyer  slant  tube 
manometer  topped  onto  the  gas  line  near  the  furnace. 

Gas  quality  was  obtained  by  gas  chromatograph  analysis  at 
the  CSU  Chemistry  Department.  One  liter  sample  bottles 
were  connected  to  the  gas  line  and  a  small  quantity  of  gas 
was  allowed  to  run  through  the  bottle  and  outside.  The 
bottles  are  double  ended  with  vacuum  tight  valves  at  each 
end.  After  purging  the  bottle,  a  vacuum  pump  is  used  for 
several  hours  to  pull  remaining  air  and  gas  out  of  the 
bottle  before  a  gas  sample  is  allowed  to  enter  the 
evacuated  bottle. 

The  sample  bottles  were  delivered  to  the  Mass  Spectroscopy 
Lab  at  the  Colorado  State  University  Chemistry  Department 
for  gas  chr  atograph  analysis  to  determine  the  constituent 
makeup  and  ■’nare  to  a  natural  gas  standard  for 
determining  ting  value.  These  samples  were  taken  an 
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average  of  every  two  days  when  many  test  runs  were  being 
made. 


Test  Procedure 

The  proposed  testing  matrix  for  determining  the  effect 


outside  combustion  air  has 
follows : 

on  furnace 

efficiency  was  as 

Outside  Air 

Temperature  Ranaes  of 

Vent  Ian 

D.uty  cycle  Percent 

-10,0,10,25,40,55 

OFF 

10,20,30,50,75,100 

LOW 

10,20,30,50,75,100 

HIGH 

10,20,30,50,75,100 

Indoor  Tempered 

AiX 

OFF 

10,20,30,50,75,100 

LOW 

10,20,30,50,75,100 

HIGH 

10,20,30,50,75,100 

Testing  was  usually  scheduled  so  that  one  complete  set  of 
duty  cycles;  i.e.,  10  through  100  percent,  for  a  given  test 
situation  could  be  completed  in  one  day.  This  was  in  fact 
rarely  possible  because  of  the  inability  to  control  the 
outside  air  temperature  over  the  length  of  this  test  run. 
The  duty  cycle  was  selected  to  be  30  minutes  as  previous 
tests  had  shown  the  furnace  operation  was  well  stabilized 
within  this  time  period.  The  percent  of  duty  cycle 
represents  the  actual  time,  as  a  percent  of  30  minutes, 
that  the  burner  was  actually  on.  The  duty  cycle  begins 
when  the  "thermostat"  activates  turning  the  burner  solenoid 
on  and  allowing  gas  to  flow  to  the  burner.  For  the  test 
series,  a  clock  timer  was  used  as  the  "thermostat"  contacts 
so  that  the  burner-on  times  could  be  easily  controlled  and 
adjusted. 

Each  30-minute  test  sequence  included  two  preconditioning 
sequences  of  the  same  burner-on  time  as  the  test  sequence, 
run  immediately  prior  to  the  test  sequence.  Testing 
usually  proceeded  beginning  with  the  10  percent  duty  cycle 
and  working  towards  the  100  percent  duty  cycle.  Again, 
because  of  the  inability  to  control  the  outside  vent  air 
temperature,  this  was  not  always  the  case.  With  the  two 
preconditioning  cycles  used,  the  test  cycle  would  reach  a 


uniform  response  regardless  of  the  test  cycle  run 
previously.  Attempts  were  made,  however,  to  not  run  the 
shorter  duty  cycles  such  as  10  and  20  percent  following  a 
100  percent  duty  cycle  run. 

During  a  test  run,  the  Doric  Data  Acquisition  System  was 
set  to  scan  the  30  data  channels  continuously  at  a  rate  of 
approximately  one  scxan  every  2.25  seconds,  resulting  in 
approximately  800  individual  data  points  for  each  data 
channel.  At  the  same  time,  all  manual  readings  were  taken 
and  the  scanivalue  stepping  switch  was  manually  stepped  to 
sample  all  air  flow  channels  several  times  during  the  test 
sequence.  The  test  was  concluded  when  the  thermostat 
switch  once  again  engaged  signalling  the  beginning  of  the 
next  30-minute  cycle.  A  known  30-minute  span  was  recorded 
since  the  thermostat  contact  could  not  provide  the  end 
signal  in  the  100  percent  case. 


Baseline  (Indoor  Tempered  Air) 

The  baseline  test  for  comparing  the  effect  of  outdoor  air 
on  furnace  efficiency  was  a  complete  series  of  the  burner 
duty  cycles  run  with  inside  or  tempered  air  as  the  source 
of  combustion  and  dilution  air  in  the  furnace  room.  The 
covers  were  placed  over  the  outside  air  vents  and  the 
covers  removed  from  the  inside  air  vents  to  accomplish 
this.  This  allowed  laboratory-tempered  air  to  enter  the 
furnace  room  —  the  same  air  source  as  is  used  for  the  room 
air  circuit.  This  series  of  tests  included  the  two-stack 
fan-induced  draft  speeds  to  simulate  moderate  and  heavy 
wind.  One  of  the  primary  requirements  for  running  this 
test  was  to  define  the  furnace  performance  for  the 
calculations  to  be  used  in  arriving  at  comparisons  of  the 
annual  fuel  utilization  for  inside  and  outside  combustion 
air . 


Outside  Combustion  Air 

The  inside  air  vents  were  blocked  and  the  outside  air  vents 
were  opened  to  allow  unheated  outside  air  to  provide  all 
combustion  and  dilution  air  to  the  furnace.  Table  A-l 
shows  that  the  complete  series  of  tests  run  for  the 
baseline  was  to  be  duplicated  for  each  indicated  outdoor 
air  temperature.  This  is  the  hoped  for  situation,  but  is 
is  recognized  that  some  portions  of  the  matrix  simply  will 
not  be  filled  in  due  to  actual  outdoor  conditions.  This  is 


59 


particularly  true  of  the  colder  temperatures  where  there 
will  be  a  limited  time  to  get  all  data  possible. 
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Simulation 

Upon  obtaining  the  furnace  performance  over  as  wide  a  range 
of  conditions  as  possible,  as  described  by  Table  A,  a 
furnace  and  nouse  model  was  formed  for  use  in  a  computer 
simulation  to  determine  annual  fuel  utilization  numbers. 

Two  house  models  were  formed,  one  with  outside  combustion 
air  to  the  furnace  and  the  other  with  tempered  indoor  air 
to  the  furnace.  A  penalty  is  included  in  the  house  model 
that  uses  indoor  combustion  air  that  takes  into  account  the 
additional  air  exfiltration  from  the  house  caused  by 
furnace  operation.  Wind  effects  are  included  in  this  model 
as  they  have  a  significant  bearing  on  the  total  stack  flow 
rate,  and  flow  rates  for  burner-on  and  burner-off  are 
included . 

A  computer  program  known  as  TRNSYS  has  been  selected  to  do 
the  performance  modeling  based  on  prior  experience  with  the 
program  at  CSU.  This  program  was  developed  by  the  Solar 
Energy  Laboratory  at  the  University  of  Wisconsin,  Madison. 
Although  it  was  originally  developed  for,  and  is  primarily 
used  on  solar  heating  and  cooling  simulations,  there  is 
nothing  inherent  within  the  program  which  would  exclude  its 
use  for  the  furnace  simulation.  In  fact,  quite  the 
opposite  is  true;  it  is  ideally  suited  for  obtaining  the 
overall  fuel  utilization  figures  for  the  furnace  models. 

Actual  meteorological  data  contained  on  "Typical 
Meteorological  Year"  ( " TMY " )  tapes  are  used  in  the 
simulation  to  control  the  house  and  furnace  model.  Three 
cities  from  among  the  listing  of  those  available  were 
selected  to  provide  a  wide  range  of  conditions  under  which 
the  house/furnace  model  would  operate.  These  cities  are 
Boston,  MA;  Great  Falls,  MO;  and  Albuquerque,  NM.  These 
three  cities  represent  a  fairly  wide  range  in  both  the 
minimum  design  temperature  and  the  degree  day  load  that  the 
furnace  must  meet. 

The  house  size  and  load  factors  were  approximately  matched 
to  the  furnace  with  a  50  percent  oversizing  of  the 
furnace.  This  has  been  a  fairly  typical  practice  in  the 
past  since  the  cost  of  additional  Btu  input  was  small  and 
natural  gas  was  cheap.  It  is  important  that  the  load 
reflect  some 
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oversizing  of  the  furnace  since  this  results  in  the 
inefficiencies  of  short  burn  cycles. 


The  computer  simulation  was  set  up  to  run  in  3-minute 
increments  for  the  months  of  October  through  April.  While 
some  small  load  may  occur  outside  these  months  in  some  of 
the  locations,  it  was  felt  that  the  additional  computer 
costs  involved  to  process  this  relatively  small  part  of  the 
load  did  not  justify  running  beyond  the  months  selected. 


APPENDIX  B 

OUTLINE  OF  KEY  DATA  FOR  EACH  CONDITION 
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Indoor 

20t 

Indoor 
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Indoor 
50 1 

Indoor 
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Indoor 

1001 

Indoor 
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Two  different  equations  must  be  used  to  approximate  the 
nonlinear  behaviorism  of  thermocouples  in  the  calculation 
of  thermopile  temperature  differentials. 

First,  the  reference  thermocouple  at  one  end  of  the 
thermopile  must  be  known  and  the  temperature  read.  This 
known  temperature  is  then  converted  to  an  equivalent 
thermocouple  EMF  relative  to  0°C  by  using  an  8th  order 
fit  equation: 

EMFref  T  =  (3.8740773840  x  10  x  T  +  3.3190198092  x  10"2  x  T2 
+  2.0714183645  x  10~4  x  T3  -  2.1945834823  x  10~6  x  T4 
+  1.1031900550  x  10~8  x  T5  -  3.0927581898  x  lO"11  x  T6 
+  4.5653337165  x  10~14  x  T7  -  2.7616878040  x  10~17  x  T8)  x  10 

where  EMF  is  in  milivolts  and  T  is  degrees  Celsius. 

The  recorded  thermopile  voltage  is  then  divided  by  the 
known  number  of  junction  pairs  making  up  the  series 
thermopile.  This  results  in  an  average  thermocouple  EMF 
difference  between  the  reference  thermocouple  and  the 
opposite  and  of  the  thermopile.  Assuming  that  the 
reference  thermocouple  temperature  is  at  the  colder  end  of 
the  thermopile  (which  is  the  way  it  was  done  in  the 
computer)  the  reference  EMF  is  added  to  the  average 
thermopile  EMF  to  obtain  an  overall  EMF  referenced  to  zero 
degrees  C. 


EMF 


EMF  +  EMF 

REF  T  AVG  T-PILE 


where  EMF 


overall 

T-PILE  Milli-Volts 


AVG  T-PILE 


N—  Junction  Pairs 


HI 


The  EMF0V^ran  is  now  converted  to  a  temperature  relative 
to  ooc  using  a  fourth  order  fit  equation: 

Toverall  =  2.5661297  x  10  x  E  -  6.1954869  x  10~1  x  e2 

+  2.2181644  x  10~2  x  E3  -  3.55000  x  10~4  x  E4 
where  E  *  EMFoverali 

T  *  Temperature  in  degrees  Celsius 
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The  thermopile  average  AT  is  now  the  difference  between 
Toverall  and  the  original  reference  temperature. 

aTt-pile  =  Toverall  ~  Treference 


Flow  rates  within  the  system  were  measured  with  pitot 
static  tubes  for  which  a  pressure  differential  recorded  as 
a  voltage.  The  voltage  was  converted  to  a  pressure 
differential  from  the  calibration  performed  on  the  Barocel. 
The  voltage-pressure  relationship  is  a  straight  line 
equation  of  the  form 

P  =  A (VOLTS)  +  B 

where  A  and  B  are  constants  from  the  calibration. 

Possible  zero  shift  within  the  Barocel  is  allowed  for  by 
measuring  one  position  in  which  a  short  circuit  tubing  is 
installed  so  that  a  known  zero  pressure  exists  across  the 
pressure  sensing  unit.  Any  zero  shift  detected  is  then 
added  or  subtracted  as  appropriate. 

The  local  velocity  can  be  calculated  when  the  pressure 
difference  on  the  pitot  static  tubes  is  known. 

VL  =  K 

and  the  local  flow  rate  is  then  the  velocity  times  the 
cross-section  area  of  the  flow: 

vL  =  ak  l/2A£_a 

'  PL 

where : 

Vl  =  local  volume  flow  rate  FT3/sec 
A  =  cross-sectional  area  of  flow  measurement  FT2 
K  =  flow  monitor  coefficient  from  calibration 
Ap  =  pitot  static  pressure  differential  lbf/ft2 
Pl  =  local  density  of  the  flow  medium  lbm/ft3 
g  =  gravity  constant  ft/sec2 


eL  is  related  to  Ps,  standard  density,  by  the  perfect 
gas  law  and  knowing  local  barometric  pressure  and  the 
temperature  at  the  pitot  static  tube,  T^: 
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psPLTs 

TLPs 


where  subscripts  s  refer  to  standard  conditions  and  L  refer 
to  local  conditions 


T  is  absolute  temperature 
P  is  barometric  pressure 


therefore: 


In  order  to  directly  add  and  subtract  flow  rates  as  measure 
by  the  two  room  air  circuit  monitors,  the  local  flow  rates 
are  converted  to  standard  conditions  by: 


v 

s 


again  using  the  relationship  for  Pl 


V 

s 


V, 


T  p 
s 

tlp 


L 

s 


or 


knowns  data 
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The  heat  flow  rate  is  calculated  .•  owing  the  temperature 
increase  from  the  earlier  T  thermopile  calculation. 

Q  =  MC  AT 
P 

where : 

Q  =  heat  flow  rate,  Btu/sec 

M  =  mass  flow  rate  Vsps,  lfc>m/sec 

Cp  =  specific  heat  of  flow  medium,  Btu/°F 

AT  =  thermopile  temperature  difference,  °F 

The  incremental  heat  flow  is: 

Ql-2  =  MCp  T  ( 1 2  -  tl) 

where  t2  -  ti  is  the  incremental  time  of  one  data  scan 
in  seconds. 


All  the  incremental  heat  flows  are  summed  in  the  computer 
for  the  test  resulting  intotal  heat  flows  for  the  test. 


The  heating  value  for  natural  gas  obtained  by  gas 
chromatograph  analysis  is  referred  to  standard  conditions. 
To  obtain  the  gas  input  Btu's  for  a  test,  this  value  must 
be  converted  to  test  conditions. 


where : 

3 

Hs  =  is  the  standard  heating  value  of  the  gas,  Btu's/Ft 
p  =  barometric  pressure,  standard  and  local 
Ts  =  standard  temperature 

Tl  =  temperature  of  the  gas  delivered  at  the  furnace 
then: 


GAS  INPUT  =  HL  .  M 


where  M  =  meter  reading,  FT3 

Barometric  pressure  was  provided  by  the  Colorado 
Climatological  Office  at  Colorado  State  University,  which 
operates  a  weather  station  at  CSU. 

Gas  Input  =  E  all  heat  flows  +  losses 

The  energy  balance  check  was  done  in  the  computer  using  the 
above  equations  and  recorded  data  to  calculate  all  heat 
flows  in  the  system. 

Jacket  and  duct  losses  in  the  system  were  not  measured  but 
were  only  estimated  to  be  1  percent  of  the  sum  of  all  other 
heat  flows.  In  steady  state  (100  percent  burner  cycle)  the 
heat  balance  was  within  1  or  2  percent  for  most  test  runs. 
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SIMULATION  OF  FURNACE  PERFORMANCE 


72 


Simulation  Program 


The  furnace  simulations  were  performed  using  a  program 
called  "TRNSYS",  a  transient  simulation  program  <?'  sloped 
by  the  University  of  Wisconsin,  Madison.  This  si  .  ’.ation 
program  was  originally  developed  for  use  in  simulating 
solar  energy  systems,  but  with  the  provision  of  a  full 
array  of  outside  heat  sources,  it  is  ideally  suited  for 
simulation  of  heating  systems  of  a  nonsolar  nature. 

"TRNSYS"  is  a  modular  program  with  subroutines  already 
developed  for  the  various  components  in  the  heating  system, 
such  as  the  load  model,  which  is  a  standard  UA  (house  load) 
model  with  single  mode  capacities.  Also  available  is  a 
thermostat  with  optional  setback  and  adjustable  hysteresis. 
A  model  for  the  furnace  was  developed  based  on  the  actual 
measured  performance  of  the  furnace  under  test.  This 
furnace  model  was  then  incorporated  into  the  "TRNSYS" 
program .. 


Load  Model 

A  standard  UA  load  model  was  used  to  predict  a  load  for 
each  selected  geographic  location.  The  UA  was  chosen  in 
each  location  so  that  the  furnace  was  about  50  percent 
oversize  for  the  load.  This  sizing  is  based  on  a  minimum 
design  temperature  for  the  particular  location,  and  the 
nameplate  output  of  the  furnace. 

When  heated  room  air  is  used  as  the  combustion  and  dilution 
air  source  for  the  furnace,  the  air  infiltration  into  the 
house  increases,  making  the  load  greater  than  for  outside 
combustion  air.  This  additional  load  is  calculated  based 
on  the  temperature  of  the  outside  air,  whether  or  not  the 
furnace  is  operating  and  the  outdoor  wind  speed  causing 
induced  draft.  This  additional  load  is  added  to  the  UA 
load  already  being  used,  as  a  penalty  for  additional  air 
intrusion.  Some  internal  heat  generation  is  also  added  to 
the  load  model  to  simulate  occupancy  and  electrical  usage. 

A  listing  of  the  system  parameters  used  to  set  up  the  model 
simulation  are  presented  in  Table  D1 . 

The  simulation  program  was  set  up  to  operate  on  three 
minute  increments.  That  is,  every  three  minutes,  the 
entire  system  is  examined  to  determine  any  changes  in 
status  or  operation.  This  coincided  with  the  minimum 
duty  cycle  measured  on  the  test  furnace  and  was  one  of  the 


key  reasons  for  selecting  a  three  minute  increment  6f  time 
step. 


The  weather  data  used  to  drive  the  model  is  contained  on 
TMY  (Typical  Meteorological  Year)  tapes  developed  by  the 
National  Oceanographic  and  Atmospheric  Administration 
(NOAA) . 


TABLE  D  1 
SYSTEM  PARAMETERS 


Location 

Minimum 
Design 
Teaip.  °F 

Load 

UA 

BtuH/°F 

Capacitance 

Btu/°F 

Internal  Heat 
Generation 

KWH /Day 

Thenaoatat 

Setpoint 

Thermostat 

Hysteresis 

Albuquerque 

14.0 

1050 

5000 

29.5 

70* 

2 

Boston 

6.0 

920 

5000 

29.5 

70* 

2 

Graat  Falls 

-20.0 

650 

5000 

29.5 

70* 

2 

•  One  option  included  a  setback  thermostat  reducing  the  setting  to  60°F 
between  the  hours  of  11  PM  and  7  AM. 
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Results 


Weekly  summaries  of  minimum  and  maximum  outdoor 
temperature,  as  well  as  overall  furnace  output,  gas  input 
and  furnace  efficiency  are  contained  in  Tables  2D  through 
13D.  These  tables  cover  the  period  1  October  through  30 
April  where  Week  1  in  the  table  refers  to  the  first  seven 
days  in  October. 

Figures  ID  through  12D  are  histogram  plots  of  furnace  run 
time  for  each  simulation.  The  height  of  the  column 
represents  the  number  of  times  during  the  simulated 
heating  season  that  the  furnace  operated  for  the  specified 
time  period.  The  furnace  had  slightly  more  operations  of 
longer  duration  when  inside  combustion  air  was  used.  With 
a  setback  thermostat,  the  furnace  shows  some  runs  of  60 
minute  duration  which  do  not  occur  when  no  setback 
thermostat  is  used.  This  is  caused  by  changing  the 
thermostat  setting  at  7  AM  from  60of  to  70OF  which 
requires  the  furnace  to  operate  for  a  longer  period  to 
recover . 

It  is  also  interesting  to  note  that  the  only  time  the 
furnace  has  any  operating  cycles  of  3  or  6  minute  durations 
is  when  a  setback  thermostat  is  used.  These  short  duty 
cycles  occur  for  all  three  locations  with  the  setback 
thermostat  in  use. 

Table  14D  and  15D  present  the  annual  fuel  data  for  the 
simulated  furnace  in  each  geographical  location  for  both 
indoor  and  outdoor  combustion  air  and  a  percent  savings 
value.  Also  presented  for  interest  are  the  results  when  a 
setback  thermostat  is  used. 

Annual  gas  savings  can  be  realized  for  each  of  the  three 
locations  simulated  when  outside  combustion  air  as  opposed 
to  inside  combustion  air  is  used  in  the  furnace  room.  The 
savings  are  greatest  for  the  most  severe  climate  of  Great 
Falls,  MT,  where  a  6  percent  annual  fuel  reduction  occurs 
and  the  least  for  the  more  moderate  climate  of  Albuquerque, 
NM,  where  the  fuel  savings  are  3.8  percent. 

Use  of  outside  combustion  air  results  in  savings  despite 
the  fact  that  the  furnace  operates  with  slightly  less 
thermal  efficiency.  Fortunately,  the  furnace  does  not 
normally  operate  for  short  10  percent  or  20  percent  duty 
cycles  (3  or  6  minute  runs)  as  seen  on  Figures  ID  through 
6D,  since  these  are  the  duty  cycles  where  the  most  severe 


drop  in  turnace  efficiency  occurs  as  outdoor  temperature 
decreases. 


With  the  load  as  modeled,  the  furnace  operates  mostly  on  12 
to  15  minute  duty  cycles,  regardless  of  the  combustion  air 
source  or  the  use  of  a  setback  thermostat. 

The  major  change  introduced  by  the  setback  thermostat  is 
the  use  of  rather  long  duration  duty  cycles  (60  minutes) 
and  also  the  addition  of  a  few  fewer  short  duration  duty 
cycles  (3  and  6  minutes) .  The  use  of  a  setback  thermostat 
still  provides  a  savings  over  no  setback  thermostat. 
However,  the  savings  are  now  greatest  for  the  more  moderate 
climate  in  Albuquerque  and  slightly  less  for  the  more 
severe  climate  at  Great  Falls. 


THERMOSTAT  SETBACK 
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FUKNACE  PERFORMANCE  summary 

LOCATION  ;  BOSTON,  MASS.  INSIDE  AIR 

THERMOSTAT  SETBACK  :  NO 
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FIGURE  Dl.  Seasonal  Histogram  of  Furnace  Operating  Cycles  for 

Albuquerque,  NM  Using  Outside  Combustion  and  Dilution  Air 
With  No  Furnace  Setback. 
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FIGURE  D2.  Seasonal  Histogram  of  Furnace  Operating  Cycles  for 

Albuqueii^ue,  NM  Using  Inside  Combustion  and  Dilution 
With  No  Furnace  Setback. 
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FIGURE  D3.  Seasonal  Histogram  of  Furnace  Operating  Cycles  for 

Boston,  MA  Using  Outside  Combustion  and  Dilution  Air 
With  No  Furnace  Setback. 
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FIGURE  D4.  Seasonal  Histogram  of  Furnace  Operating  Cycles  for 
Boston,  MA  Using  Inside  Combustion  and  Dilution  Air 
With  No  Furnace  Setback. 
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FIGURE  D5.  Seasonal  Histogram  of  Furnace  Operating  Cycles  for 

Great  Falls,  KT  Using  Outside  Combustion  and  Dilution  Air 
With  No  Furnace  Setback. 
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FIGURE  D6.  Seasonal  Histogram  of  Furnace  Operating  Cycles  for 

Great  Falls ,  MT  Using  Inside  Combustion  and  Dilution  Air 
With  No  Furnace  Setback. 
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FIGURE  D7 .  Seasonal  Histogram  of  Furnace  Operating  Cycles  for 

Albuquerque,  NM  Using  Outside  Combustion  and  Dilution  Air 
With  A  Furnace  Setback. 
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FIGURE  D8.  Seasonal  Histogram  of  Furnace  Operating  Cycles  for 

Albuquerque,  NM  Using  Inside  Combustion  and  Dilution  Air 
With  A  Furnace  Setback. 
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FIGURE  D9.  Seasonal  Histogram  of  Furnace  Operating  Cycles  for 

Boston,  MA  Using  Outside  Combustion  and  Dilution  Air 
With  A  Furnace  Setback. 
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FIGURE  DIO.  Seasonal  Histogram  of  Furnace  Operating  Cycles  for 
Boston,  MA  Using  Inside  Combustion  and  Dilution  Air 
With  A  Furnace  Setback. 
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FIGURE  Dll.  Seasonal  Histogram  of  Furnace  Operating  Cycles  for 

Great  Falls,  MT  Using  Outside  Combustion  and  Dilution  Air 
With  A  Furnace  Setback. 


100 


NO.  OF  TIMES  SYSTEM  OPERATED 


4000 


3000 


2000 


1000 


0 


FIGURE  DL2. 


0> 

sr 

OJ 


Seasonal  Histogram  of  Furnace  Operating  Cycles  for 
Great  Falls,  MI  Using  Inside  Combustion  and  Dilution  Air 
With  A  Furnace  Setback. 
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TABLE  PI 4 

ANNUAL  FUEL  US ACE  AND  EFFICIENCY 


Inside  Coabuation  Air 

Outsl 

de  Coabust 

Ion  Air 

Theraostat 

Caa  Input 

Boos  Air 
Delivery 

Efficiency 

Css  Input 

Boon  Air 
Delivery 

Efficiency 

Location 

Type 

BTUnlO'6 

lTUiilO'* 

X 

BTUalO"6 

BTUalO'4 

Z 

- - j- 

Albuquerqoe , 

m 

No  Night 
Setback 

187.5 

131.4 

70.1 

180.3 

124.0 

68.8 

Night 

Setback 

161.6 

114.0 

70.3 

155.7 

107.8 

69.2 

i 

w 

ii 

No  Night 
Setback 

192.6 

134.3 

69.7 

183.2 

125.0 

68.2 

Night 

Setback 

169.6 

119.2 

70.3 

161.4 

111.1 

68.8 

Great  Falls, 

m 

No  Night 
Setback 

163.8 

113.7 

69.4 

153.9 

103.1 

67.0 

Night 

Setback 

146.6 

103.0 

70.3 

137.9 

93.5 

67.8 

TABLE  D15 

COMPUTED  SAVINGS  FOB  THE  VARIOUS  TEST  CASES 


Saving.  1 

Load  Design  Paraneters 

Location 

A 

B 

C 

D 

Design 

Teaperature 

*F 

Actual 

T  Kin 
•F 

Load 

UA 

BTUH/#I 

Albuquerque 

3.84 

3.65 

13.81 

13.84 

14.0 

8.9 

1050 

Boston 

4.88 

4.83 

11.94 

11.90 

6.0 

6.0 

920 

Creat  Palls 

6.04 

5.93 

10.50 

10.40 

-20.0 

-17.0 

650 

, Inside-Out  aide  v 

A.  Inside  vs.  outside  coabustloo  sir -no  setback  thensostst  ( - \nti6e - * 

b.  Inside  vs.  outside  coabustlon  air-vith  setback  thermostat  (  "  ) 

C.  Setback  vs.  no  setback  thermostat-inside  air 

D.  Setback  vs.  no  setback  theraostat -out  aide  air 


.No  Setback-Setback. 
1  No  Setback  ; 

(  "  ) 
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US  Military  Academy 
ATTN:  Dept  of  Mechanics 
West  Point,  NY  10996 

US  Military  Academy 
ATTN:  Library 
West  Point,  NY  10996 


Comnander,  TRADOC 
Office  of  the  Engineer 
ATTN:  ATEN-FE-U 
Ft  Monroe,  VA  23651 

A F  Civil  Engr  Center/XRL 
Tyndall  AFB,  FL  32401 


HQDA  (DALO-TSE-F) 

WASH  DC  20314 

HQDA  (DAEN-ASI-L)  (2) 
WASH  DC  20314 

HQDA  (DAEN-MPO-B) 

WASH  DC  20314 

HQDA  (DAEN-MPR-A) 

WASH  DC  20314 

HQDA  (DAEN-MPO-U) 

WASH  DC  2031 4 


Naval  Facilities  Engr  Command 
ATTN:  Code  04 
200  Stovall  St. 

Alexandria,  VA  22332 

Defense  Documentation  Center 
ATTN:  TCA  (12) 

Cameron  Station 
Alexandria,  VA  223  1  4 

Commander  and  Director 
USA  Cold  Regions  Research  Engineering 
Laboratory 
Hanover,  NH  037  55 


HQDA  (DAEN-MPZ-A) 
WASH  DC  20314 

HQDA  (DAEN-MPZ-E) 
WASH  DC  20314 

HQDA  (DAEN-MPZ-G) 
WASH  DC  20314 

HQDA  (DAEN-RDM) 

WASH  DC  20314 

HQDA  (DAEN-RDL) 

WASH  DC  20314 

Director,  USA-WES 
ATTN:  Library 
P.0.  Box  631 
Vicksburg,  MS  39181 


FORSCOM 
ATTN:  AFEN 

Ft  McPherson,  GA  30330 
FORSCOM 

ATTN:  AFEN-FE 
Ft  McPherson,  GA  30330 

Offlcer-in-Charge 
Civil  Engineering  Laboratory 
Naval  Construction  Battalion  Center 
ATTN:  Library  (Code  L08A) 

Port  Hueneme,  CA  93043 

Commander  and  Director 
USA  Construction  Engineering 
Research  Laboratory 
P.0.  Box  4005 
Champaign,  IL  61820 


Commander,  TRADOC 
Office  of  the  Engineer 
ATTN:  ATEN 
Ft.  Monroe,  VA  23651 


Commanding  General,  3d  USA 

ATTN:  Engineer 

Ft.  McPherson,  GA  303  3  0 
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Commanding  General,  5th  USA 

ATTN:  Engineer 

Ft  Sam  Houston,  TX  78234 

AFCE  Center 
Tyndall  AFB,  FL  32403 

Commander,  0 ARC  CM 
Director,  Installation 
and  Services 
5001  Eisenhower  Ave. 

Alexandria,  VA  22333 

Commander,  DARCGM 

ATTN:  Chief,  Engineering  Div. 

5001  Eisenhower  Ave 
Alexandria,  VA  22333 

Air  Force  Weapons  Lab/AFWL/DE 
Chief,  Civil  Engineering 
Research  Division 
Kirtland  AFB,  NM  87117 

Strategic  Air  Command 
ATTN:  DSC/CE  (DEEE) 

Offutt  AFB,  NE  68112 

Headquarters  USAF 

Directorate  of  Civil  Engineering 

AF/PREES 

Bolling  AFB,  Washington,  DC  20333 

Strategic  Air  Command 
Engineeri ng 
ATTN:  Ed  Morgan 
Offutt  AFB,  NE  68113 

USAF  Institute  of  Technology 
AF  IT/D  ED 

Wright  Patterson  AFB,  OH  45433 

Air  Force  Weapons  Lab 
Technical  Library  (DOUL) 

Kirtland  AFB,  NM  87117 


Chief,  Naval  Facilities 
Engineer  Command 
ATTN:  Chief  Engineer 
Department  of  the  Navy 
Washington,  DC  20350 

Commander 

Naval  Facilities  Engineering  Qnd 
200  Stovall  St 
Alexandria,  VA  22332 

Commander 

Naval  Facilities  Engr  Qnd 
Western  Division 
Box  727 

San  Bruno,  CA  94066 

Civil  Engineering  Center 
ATTN:  Moreell  Library 
Port  Hueneme,  CA  93043 

Commandant  of  the  Marine  Corps 
tfQ,  US  Marine  Corps 
Washington,  DC  20380 

National  Bureau  of  Standards  (4) 
Materials  &  Composites  Section 
Center  for  Building  Technology 
Washington,  DC  20234 

Assistant  Chief  of  Engineer 
Rn  IE  668,  Pentagon 
Washington,  DC  20310 

The  Army  Library  ’(ANRAL-R) 

ATTN:  Army  Studies  Section 
Room  1A  518,  The  Pentagon 
Washington,  DC  2031  0 

Comma  nder- In-Chief 

USA,  Europe 

ATTN:  AEAEN 

APO  New  York,  NY  09403 
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Commander 

USA  Foreign  Science  and 
Technology  Center 
220  8th  St.  N.E. 
Charlottesville,  VA  22901 

Commander 

USA  Science  &  Technology 
Information  Team,  Europe 
APO  New  York,  NY  09710 

Commander 

USA  Science  &  Technology 
Center  -  Far  East  Office 
APO  San  Francisco,  CA  96328 

Commanding  General 

USA  Engineer  Comnand,  Europe 

APO  New  York,  NY  09403 

Deputy  Chief  of  Staff 
for  Logistics 
US  Army,  The  Pentagon 
Washington,  DC  20310 

Commander,  TRAD0C 
Office  of  the  Engineer 
ATTN:  Chief,  Facilities 
Engineering  Division 
Ft  Monroe,  VA  23651 

Commanding  General 
USA  Forces  Command 
Office  of  the  Engineer 
(AFEN-FES) 

Ft  McPherson,  GA  30330 

Commanding  General 
USA  Forces  Command 
ATTN:  Chief,  Facilities 
Engineering  Division 
Ft  McPherson,  GA  3033  0 


Commander 

USA  Support  Command,  Hawaii 
Fort  Shatter,  HI  96858 

Commander 

Eighth  US  Army 

APO  San  Francisco  96301 

Commander 

US  Army  Facility  Engineer 
Activity  -  Korea 
APO  San  Francisco  96301 

Commander 

US  Army,  Japan 

APO  San  Francisco,  CA  96343 

Facilities  Engineer 

Fort  Bel  voir 

Fort  Belvoir,  VA  22060 

Facilities  Engineer 

Fori  Benning 

Fort  Benning,  GA  31905 

Facilities  Engineer 
Fort  Bliss 

Fort  3l1ss,  TX  79916 

Facilities  Engineer 
Carlisle  Barracks 
Carlisle  Barracks,  PA  17013 

Facilities  Engineer 

Fort  Chaffee 

Fort  Chaffee,  AR  .  72902 

Facilities  Engineer 
Fort  Dix 

Fort  Dix,  NJ  08640 

Facilities  Engineer 
Fort  Eustis 

Fort  Eustis,  VA  23604 


Commanding  General,  1st  USA 

ATTN:  Engineer 

Ft  George  G.  Meade,  MD  20755 


Facilities  Engineer 
Fort  Gordon 

Fort  Gordon,  GA  30905 


Facilities  Engineer 
Fort  Story 

Fort  Story,  VA  23459 


Facilities  Engineer 

Fort  Hamilton 

Fort  Hamilton,  NY  11252 

Facilities  Engineer 

Fort  A  P  Hill 

Bowling  Green,  VA  22427 

Facilities  Engineer 

Fort  Jackson 

Fort  Jackson,  SC  29207 

Facilities  Engineer 
Fort  Knox 

Fort  Knox,  KY  40121 

Facilities  Engineer 
Fort  Lee 

Fort  Lee,  VA  23801 

Facilities  Engineer 

Fort  McClellan 

Fort  McClellan,  AL  36201 

Facilities  Engineer 

Fort  Monroe 

Fort  Monroe,  VA  23651 

Facilities  Engineer 

Presidio  of  Monterey 

Presidio  of  Monterey,  CA  93940 

Facilities  Engineer 
Fort  Pickett 
Black stone,  VA  23824 

Facilities  Engineer 

Fort  Rucker 

Fort  Rucker,  AL  36362 


Facilities  Engineer 
Kansas  Army  Ammunition  Plant 
Parsons,  KS  67357 

Facilities  Engineer 

Lone  Star  Amy  Ammunition  Plant 

Texarkana,  TX  75501 

Facilities  Engineer 
Plcatlnny  Arsenal 
Dower,  NJ  07801 

Facilities  Engineer 
Louisiana  Army  Ammunition  Plant 
Shreveport,  LA  71130 

Facilities  Engineer 
Milan  Army  Ammunition  Plant 
Milan,  TN  38358 

Facilities  Engineer 
Pine  Bluff  Arsenal 
Pine  Bluff,  AR  71601 

Facilities  Engineer 
Radford  Army  Ammunition  Plant 
Radford,  VA  24141 

Facilities  Engineer 
Rock  Island  Arsenal 
Rock  Island,  IL  61201 

Facilities  Engineer 
Rocky  Mountain  Arsenal 
Denver,  CO  80340 

Facilities  Engineer 
Scranton  Army  Ammunition  Plant 
156  Cedar  Avenue 
Scranton,  PA  18503 


Facilities  Engineer 
Fort  Sill 

Fort  Sill,  OK  73503 
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Facilities  Engineer 
Tobyhanna  Army  Depot 
Tobyhanna,  PA  18466 


Facilities  Engineer 
Tooele  Army  Oepot 
Tooele,  UT  84074 

Facilities  Engineer 
Arlington  Hall  Station 
400  Arlington  Blvd 
Arlington,  VA  22212 

Facilities  Engineer 
Cameron  Station,  Bldg  17 
5010  Duke  Street 
Alexandria,  VA  22314 

Facilities  Engineer 

Sunny  Point  Military  Ocean  Terminal 

Southport,  NC  28461 

Facilities  Engineer 
US  Military  Academy 
West  Point  Reservation 
West  Point,  NY  10996 

Facilities  Engineer 

Fort  Ritchie 

Fort  Ritchie,  MO  21719 

Facilities  Engineer 
Army  Materials  &  Mechanics 
Research  Center 
Watertown,  MA  02172 

Facilities  Engineer 
Ballistics  Missile  Advanced 
Technology  Center 
P.0.  Box  1500 
Huntsville,  AL  35807 

Facilities  Engineer 
Fort  Walnwrlght 
172d  Infantry  Brigade 
Fort  Walnwrlght,  AK  99703 

Facilities  Engineer 

Fort  Greely 

Fort  Greely,  AK  98733 

Facilities  Engineer 

Fort  Richardson 

Fort  Richardson,  AK  99605 


Facilities  Engineer 
Harry  Diamond  Laboratories 
2800  Powder  Mill  Rd 
Adel  phi,  MD  20783 

Facilities  Engineer 
Fort  Missoula 
Missoula,  MT  59801 

Facilities  Engineer 
New  Cumberland  Army  Depot 
New  Cumberland,  PA  17070 


Facilities  Engineer 
Oakland  Army  Base 
Oakland,  CA  94626 

Facilities  Engineer 
Vint  Hill  Farms  Station 
Warrentown,  VA  22186 

Facilities  Engineer 

Twin  Cities  Army  Ammunition  Plant 

New  Brighton,  MN  55112 

Facilities  Engineer 
Volunteer  Army  Ammunition  Plant 
Chattanooga ,  TN  37401 

Facilities  Engineer 
Watervllet  Arsenal 
Watervllet,  NY  12189 

Facilities  Engineer 
St  Louis  Area  Support  Center 
Granite  City,  IL  62040 

Facilities  Engineer 

Fort  Monmouth 

Fort  Moranouth,  NJ  07700 

Facilities  Engineer 
Redstone  Arsenal 
Redstone  Arsenal,  AL  35809 
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Facilities  Engineer 
Detroit  Arsenal 
Warren,  MI  48039 

Facilities  Engineer 

Aberdeen  Proving  Ground 

Aberdeen  Proving  Ground,  MD  21005 

Facilities  Engineer 
Jefferson  Proving  Ground 
Madison,  IN  47250 

Facilities  Engineer 
Dugway  Proving  Ground 
Dugway,  UT  84022 

Facilities  Engineer 
Fort  McCoy 
Sparta,  UI  54656 

Facilities  Engineer 

White  Sands  Missile  Range 

White  Sands  Missile  Range,  fW  88002 

Facilities  Engineer 
tana  Proving  Ground 
tana,  AZ  85364 

Facilities  Engineer 
Natick  Research  &  Dev  Ctr 
Kansas  St. 

Natick,  MA  01760 

Facilities  Engineer 
Fort  Bragg 

Fort  Bragg,  NC  28307 

Facilities  Engineer 

Fort  Campbell 

Fort  Campbell,  KY  42223 

Facilities  Engineer 

Fort  Carson 

Fort  Carson,  CO  80913 

Facilities  Engineer 
Fort  Drum 

Watertown,  NY  13601 


Facilities  Engineer 
Fort  Hood 

Fort  Hood,  TX  76544 

Facilities  Engineer 
Fort  Indiantown  Gap 
Annvllle,  PA  17003 

Facilities  Engineer 
Fort  Lewis 

Fort  Lewis,  WA  98433 

Facilities  Engineer 

Fort  MacArthur 

Fort  MacArthur,  CA  90731 

Facilities  Engineer 

Fort  McPherson 

Fort  McPherson,  GA  3  033  0 

Facilities  Engineer 

Fort  George  G.  Meade 

Fort  George  G.  Meade,  MD  20755 

Facilities  Engineer 
Fort  Polk 

Fort  Polk,  LA  71459 

Facilities  Engineer 
Fort  Riley 

Fort  Riley,  KS  66442 

Facilities  Engineer 

Fort  Stewart 

Fort  Stewart,  GA  31312 

Facilities  Engineer 
Indiana  Army  Ammunition  Plant 
Charlestown,  IN  47111 

Facilities  Engineer 
Joliet  Army  Ammunition  Plant 
Joliet,  IL  60436 

Facilities  Engineer 
Anniston  Army  Depot 
Anniston,  AL  36201 
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Facilities  Engineer 
Corpus  Christi  Army  Depot 
Corpus  Christi,  TX  78419 

Facilities  Engineer 
Red  River  Army  Depot 
Texarkana,  TX  75501 

Facilities  Engineer 
Sacramento  Army  Depot 
Sacramento,  CA  95813 

Facilities  Engineer 
Sharpe  Army  Depot 
Lathrop,  CA  953  3  0 


Facilities  Engineer 

Gulf  Output 

New  Orleans,  LA  70146 

Facilities  Engineer 

Fort  Huachuca 

Fort  Huachuca,  AZ  86513 

Facilities  Engineer 
Letterkenny  Army  Depot 
Chambersburg,  PA  17201 

Facilities  Engineer 
Michigan  Army  Missile  Plant 
Warren,  MI  48089 


Facilities  Engineer 
Seneca  Army  Depot 
Romulus,  NY  14541 

Facilities  Engineer 
Fort  Ord 

Fort  Ord,  CA  93941 

Facilities  Engineer 

Presidio  of  San  Francisco 

Presidio  of  San  Francisco,  CA  94129 

Facilities  Engineer 

Fort  Sheridan 

Fort  Sheridan,  IL  60037 

Facilities  Engineer 
Holston  Army  Ammunition  Plant 
Kingsport,  TN  37662 


COL  E.C.  Lussier 
FI tz Simons  Army  Med  Center 
ATTN:  HSF-DFE 
Denver,  CO  80240 

US  Army  Engr  Dist,  New  York 
ATTN:  NANEN-E 
26  Federal  Plaza 
New  York,  NY  10007 

USA  Engr  Dist,  Baltimore 
ATTN:  Chief,  Engr  Div 
P.0.  Box  1715 
Baltimore,  MD  21203 

USA  Engr  Dist,  Charleston 
ATTN:  Chief,  Engr  Div 
P.0.  Box  919 
Charleston,  SC  29402 


Facilities  Engineer 
Baltimore  Output 
Baltimore,  MD  21222 


USA  Engr  Dist,  Detroit 
P.0.  Box  1027 
Detroit,  MI  48231 


Facilities  Engineer 

Bayonne  Military  Ocean  Terminal 

Bayonne,  NJ  07002 

Facilities  Engineer 

Bay  Area  Military  Ocean  Terminal 

Oakland,  CA  94626 


USA  Engr  Dist,  Kansas  City 
ATTN:  Chief,  Engr  Div 
700  Federal  Office  Bldg. 
601  E.  12th  St 
Kansas  City,  MO  64106 
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USA  Engr  Dist,  Qnaha 
ATTN:  Chief,  Engr  Div 
7410  USOP  and  Courthouse 
215  N.  17th  St 
Onaha,  NE  68102 

USA  Engr  Dist,  Fort  Worth 
ATTN:  Chief,  SWF ED- D 
P.0.  Box  17300 
Fort  Worth,  TX  76102 

USA  Engr  Dist,  Sacramento 
ATTN:  Chief,  SPKED-D 
650  Capitol  Mall 
Sacramento,  CA  95814 

USA  Engr  Dist,  Far  East 
ATTN:  Chief,  Engr  Div 
APO  San  Francisco,  CA  96301 


USA  Engr  Dist,  Japan 

APO  San  Francisco,  CA  96343 

USA  Engr  Div,  Europe 
European  Div,  Corps  of  Engineers 
APO  New  York,  NY  09757 

USA  Engr  Div,  North  Atlantic 
ATTN:  Chief,  NADEN-T 
90  Church  St. 

New  York,  NY  10007 

USA  Engr  Div,  South  Atlantic 
ATTN:  Chief,  SAEN-TE 
510  Title  Bldg 
30  Pryor  St,  SW 
Atlanta,  GA  30303 

USA  Engr  Dist,  Mobile 
ATTN:  Chief,  SAMEN-C 
P.0.  Box  2288 
Mobile,  AL  36601 

USA  Engr  Dist,  Louisville 
ATTN:  Chief,  Engr  Div 
P.0.  Box  59 
Louisville,  KY  40201 


USA  Engr  Div,  Norfolk 
ATTN:  Chief,  NAOEN-D 
803  Front  Street 
Norfolk,  VA  23510 

USA  Engr  Div,  Missouri  River 
ATTN:  Chief,  Engr  Div 
P.0.  Box  103  Downtown  Station 
Omaha,  NE  68101 

USA  Engr  Div,  South  Pacific 
ATTN:  Chief,  SPDED-TG 
630  Sansome  St,  ftn  1216 
San  Francisco,  CA  94111 

USA  Engr  Div,  Huntsville 
ATTN:  Chief,  HNDED-ME 
P.0.  Box  1600  West  Station 
Huntsville,  AL  35807 

USA  Engr  Div,  Ohio  River 
ATTN:  Chief,  Engr  Div 
P.0.  Box  1159 
Cincinnati,  Ohio  45201 

USA  Engr  Div,  North  Central 
ATTN:  Chief,  Engr  Div 
536  S.  Clark  St. 

Chicago,  IL  60605 

USA  Engr  Div,  Southwestern 
ATTN:  Chief,  SWDED-TM 
Main  Tower  Bldg,  1200  Main  St 
Dallas,  TX  75202 

USA  Engr  Dist,  Savannah 
ATTN:  Chief,  SASAS-L 
P.0.  Box  889 
Savannah,  GA  31402 

Commander 

US  Army  Facilities  Engineering 
Support  Agency 
Support  Detachment  II 
Fort  Gillem,  GA  30050 
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Commander 

US  Army  Facilities  Engr  Spt  Agency 

ATTN:  MAJ.  Brisbine 

Support  Detachment  III 

P.0.  Box  6550 

Fort  Bliss,  TX  79916 

NCOIC 

US  Army  Facilities  Engr  Spt  Agency 

Support  Detachment  III 

ATTN:  FESA-III-SI 

P.0.  Box  3  031 

Fort  Sill,  OK  73503 

NCOIC 

US  Army  Facilities  Engr  Spt  Agency 
Support  Detachment  III 
ATTN:  FESA-I II-PR 
P.0.  Box  29704 

Presidio  of  San  Francisco,  CA  94129 
NCOIC 

US  Army  Facilities  Engr  Spt  Agency 
ATTN:  FESA-III-CA 
°ost  Locator 
ort  Carson,  CO  80913 

Commander/CPT  Ryan 

US  Army  Facilities  Engr  Spt  Agency 

Support  Detachment  IV 

P.0.  Box  3  00 

Fort  Monmouth,  NJ  07703 

NCOIC 

US  Army  Facilities  Engr  Spt  Agency 

ATTN:  FESA-IV-MU 

P.0.  Box  300 

Fort  Monmouth,  NJ  07703 

NCOIC 

US  Army  Facilities  Engr  Spt  Agency 
Support  Detachment  IV 
ATTN:  FESA-IV-ST 
Stewart  Army  Subpost 
Newburgh,  New  York  12250 

NCOIC 

US  Army  Facilities  Engineering 
Support  Agency 
Support  Detachment  II 
ATTN:  FESA-I I-JA 
Fort  Jackson,  SC  29207 


NCOIC 

US  Army  Facilities  Engr  Spt  Agency 

Support  Detachment  II 

ATTN:  FESA-II-BE 

P.0.  Box  2207 

Fort  Benning  GA  31905 

NCOIC 

US  Army  Facilities  Engr  Spt  Agency 
Support  Detachment  II 
ATTN:  FESA-I I- KN 
Fort  Knox,  KY  40121 

Naval  Facilities  Engineering  Qnd 
Energy  Programs  Branch,  Code  1023 
Hoffmann  Bldg.  2,  (Mr.  John  Hughes) 
Stovall  Street 
Alexandria,  VA  22332 

Commander 

US  Army  Facilities  Engineering 
Support  Agency 
FE  Support  Detachment  I 
APO  New  York,  NY  09081 

Navy  Energy  Office 
ATTN:  W.R.  Mitchum 
Washington  DC  20350 

David  C.  Hall 

Energy  Projects  Officer 

Dept,  of  the  Air  Force 

Sacramento  Air  Logistics  Center  (AFLC) 

2852  ABG/DEE 

McClellan,  CA  9$652 

USA  Engineer  District,  Chicago 
219  S.  Dearborn  Street 
ATTN:  District  Engineer 
Chicago,  IL  60604 

Directorate  of  Facilities  Engineer 
Energy  Environmental  &  Self  Help  Center 
Fort  Campbell,  KY  42223 

Commander  and  Director 
Construction  Engineering  Research 
Laboratory 
ATTN:  COL  Clrceo 
P.0.  Box  4005 
Champaign,  IL  61820 
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Mr.  Ray  Heller 
Engineering  Services  Branch 
DFAE,  Bldg,  1950 
Fort  Sill,  OK  73503 

Commander-in-Chief 
HQ,  USAEUR 
ATTN:  AEAEN-EH-U 
APO  New  York  09403 

HQ  AFESC/RDVA 
ATTN:  Mr.  Hathaway 
Tyndall  AFB,  FL  32403 

Commander  and  Director 

Construction  Engineering  Research  Lab 

ATTN:  Library 

P.0.  Box  4005 

Champaign,  IL  61820 

HQ,  5th  Signal  Command 
Office  of  the  Engineer 
APO  New  York  09056 

HQ,  Us  Military  Community  Activity, 
Heilbronn 

Director  of  Engineering  &  Housing 
ATTN:  Rodger  D.  Romans 
APO  New  York  09176 

Commanding  General 
HQ  USATC  and  Fort  Leonard  Wood 
ATTN:  Facility  Engineer 
Fort  Leonard  Wood,  MO  65473 

SSG  Ruiz  Burgos  Andres 

O. F.E.,  HHC  HQ  Qnd  193d  Inf 
BDE 

Ft.  Clayton,  C/Z 

Energy/Environmental  Office 
ATTN:  David  R.  Nichols 
USMCA-NBG  (DEH) 

APO  New  York  09696 

Commander 

535th  Engineer  Detachment 

P. 0.  Box  300 

Fbrt  Monmouth,  NJ  077  03 


NCOIC 

535th  Engineer  Detachment,  Team  A 

ATTN:  SFC  Prenger 

P.0.  Box  224 

Fort  Knox,  KY  40121 

NCOIC 

535th  Engineer  Detachment,  Team  B 
ATTN:  SP6  Cathers 
P.0.  Box  300 

Fort  Monmouth,  NJ  07703 
NCOIC 

535th  Engineer  Detachment,  Team  C 

ATTN:  SFC  Jackson 

P.0.  Box  4301 

Fort  Eustls,  V A  23604 

NCOIC 

535th  Engineer  Detachment,  Team  D 
ATTN:  SFC  Hughes 
Stewart  Army  Subpost 
Newburg,  New  York  12550 

Commander 

Persidio  of  San  Francisco, 
California 

ATTN:  AFZM-DI/Mr.  Prugh 
San  Francisco,  CA  94129 

Facilities  Engineer 
Corpus  Christi  Army  Depot 
ATTN:  Mr.  Joseph  Canpu/Stop  24 
Corpus  Christi,  TX  78419 

Walter  Reed  Army  Medical  Center 
ATTN:  KHSWS-E/ James  Prince 
6825  16th  St.,  NW 
Washington,  DC  20012 

Commanding  Officer 

Installations  and  Services  Activity 

ATTN:  DRCIS-R I- I B 

Rock  Island  Arsenal 

Rock  Island,  IL  61299 
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Commanding  Officer 
Northern  Division  Naval 

Facilities  Engineering  Command 
Code  102  (Mr.  E.F.  Numm) 

Naval  Base 

Philadelphia,  PA  19112 

Commander,  US  Army  Facilities  Engineering 
Support  Agency 
Support  Detachment  I 
APO  New  York  09081 

HQ,  USA  Health  Services  Qnd 
Bldg.  2792 
ATTN:  HSLO-F 

Fort  Sam  Houston,  TX  78234 
HQDA 

(DAEN-MPE-E) 

WASH  DC  20314 

Commanding  Officer 
Northern  Division  Naval 

Facilities  Engineering  Conmand 
Code  10 

Naval  Base,  Building  77 
Philadelphia,  PA  19112 

Facilities  Engineer 

Fort  Leavenworth 

Fort  Leavenworth,  KS  66027 

Facilities  Engineer 

Fort  Benjamin  Harrison 

Fort  Benjamin  Harrison,  IN  46216 

Office  of  the  A&E 
ATTN:  MAJ  Johnson 
Camp  Ripley 

Little  Falls,  MN  56345 

Commander 
US  Army  Garrison 
ATTN:  HSD-FE 
Fort  Detrick,  MD  21701 


Mr.  David  White 
Defense  Audit  Service 
888  North  Sepulveda  Blvd. 

Suite  610 

El  Sequndo,  CA  90245 

Facilities  Engineer 
Bldg.  308 

Fort  Myer,  VA  22211 
NAVFAC 

ATTN:  John  Zekan 
Code  0833  Hoffmann  Building 
200  Stovall  Street 
Alexandria,  VA  22332 

HQ,  USASCH 

Director  Engineering  &  Housing 
Fort  Shafter,  HI  96858 

HQ,  WESTCOM 

ATTN:  APEN-CE 

Fort  Shafter,  HI  96858 

Headquarters  US  Army  Materiel  Development 
&  Readiness  Command 
ATTN:  Energy  Office,  DRC1S-C 
Alexandria,  VA  22333 

One  Stop  Coordinator 
Army  Corps  of  Engineers 
ATTN:  ORNED-D  (Connie  Flatt) 

P.0.  Box  1070 
Nashville,  TN  37202 

Solar  Energy  Research  Institute 
1617  Cole  Boulevard 
Golden,  CO  80401 

American  Telephone  &  Telegraph  Co. 

ATTN:  Kenneth  T.  Risberg 
222  Mt.  Airy  Road,  ftn  192B5 
Basking  Ridge,  NJ  07920 

LCDR  D.  J.  Clark 
Navy  Material  Command 
COde  08 E 

Washington,  OC  20360 


AF ESC/DEB 

ATTN:  Mr.  Fred  Beason 
Tyndall  AFB,  FL  32403 
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Office  of  Secretary  of  Defense 
Installations  &  Housing 
ATTN:  Mr.  Millard  Carr 
WASH  DC  20301 

Commandant  (G-ECV-2/65) 

ATTN:  LTC  Peck 
US  Coast  Guard  HQTRS 
2100  2nd  St.  SW 
WASH  DC  20593 

HQ  AFESC/DEB 

ATTN:  COL.  William  R.  Gaddi e 
Tyndall  AFB,  FL  32403 
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